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ABSTRACT 

We present 0.5-2 keV, 2-8 keV, 4-8 keV, and 0.5-8 keV (hereafter, soft, hard, ultra-hard, and full bands, 
respectively) cumulative and differential number counts (log A^-log S) measurements for the recently completed 
«4 Ms Chandra Deep Field-South (CDF-S) survey, the deepest X-ray survey to date. We implement a new 
Bayesian approach, which allows reliable calculation of number counts down to flux limits that are factors of 
ssL9^.3 times fainter than the previously deepest number-counts investigations. In the soft band, the most 
sensitive bandpass in our analysis, the ?s4 Ms CDF-S reaches a maximum source density of w27,800 deg"^. 
By virtue of the exquisite X-ray and multiwavelength data available in the CDF-S, we are able to measure 
the number counts from a variety of source populations (active galactic nuclei [AGNs], normal galaxies, and 
Galactic stars) and subpopulations (as a function of redshift, AGN absorption, luminosity, and galaxy morphol- 
ogy), and test models that describe their evolution. We find that AGNs still dominate the X-ray number counts 
down to the faintest flux levels for all bands and reach a limiting soft-band source density of « 14,900 deg"^, 
the highest reliable AGN source density measured at any wavelength. We find that the normal-galaxy counts 
rise rapidly near the flux limits, and at the limiting soft-band flux, reach source densities of « 12,700 deg"^ 
and make up 46 ± 5% of the total number counts. The rapid rise of the galaxy counts toward faint fluxes, 
and significant normal-galaxy contributions to the overall number counts, indicate that normal galaxies will 
overtake AGNs just below the «4 Ms soft-band flux limit and will provide a numerically significant new X-ray 
source population in future surveys that reach below the w4 Ms sensitivity limit. We show that a future ~ 10 Ms 
CDF-S would allow for a significant increase in X-ray detected sources, with many of the new sources being 
cosmologically distant (z ^ 0.6) normal galaxies. 

Subject headings: cosmology: observations — galaxies: active — galaxies: starburst — X-rays: galaxies 



1. INTRODUCTION 

Deep extragalactic X-ray surveys conducted with Chan- 
dra and XMM-Newton have resolved the vast majority of the 
0.5-10 keV cosmic X-ray background (CXRB) and have pro- 
vided substantial new insight into the X-ray point sources de- 
tected (see, e.g., Brandt & Hasinger 2005; Brandt & Alexan- 
der 2010). A fundamental quantity used to characterize the 
extragalactic X-ray source population is the cumulative X-ray 
number counts, which quantifies how the cumulative number 
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of X-ray sources per unit area, A^, increases with decreasing 
flux, S (also referred to as logA^-logS'; see, e.g., Brandt etal. 
2001; Rosati etal. 2002; Moretti etal. 2003; Bauer etal. 2004; 
Kim et al. 2007; Georgakakis et al. 2008). The number counts 
provide an observational constraint that must be considered 
when constructing physical models that describe the X-ray 
evolution of extragalactic sources in the Universe. For exam- 
ple, successful models including the supermassive black hole 
(SMBH) accretion history of the Universe (e.g., Gilli etal. 
2007; Treister et al. 2009) and the X-ray evolution of normal 
galaxies (primarily driven by X-ray binaries and hot gas; see, 
e.g., Ranalli et al. 2003, 2005) must predict X-ray source den- 
sities consistent with the observed number counts. 

Due to its high angular resolution and low background, 
Chandra is the only current X-ray observatory capable of pro- 
viding new constraining data in the ultra-deep regime (be- 
low 0.5-8 keV fluxes of a few x 10"'^ ergs cm"^ s~'), since 
the deepest Chandra surveys have now greatly surpassed the 
XMM-Newton confusion limit (see, e.g., Brandt etal. 2001; 
Giacconi etal. 2001, 2002; Alexander etal. 2003; Luo etal. 
2008; Xue et al. 201 1). As new Chandra surveys continue to 
probe the extragalactic X-ray universe to fainter depths, the 
number counts continuously rise as fainter X-ray populations 
are revealed. Naturally, it is at the faintest flux levels where 
ultradeep Chandra surveys are probing new regions of discov- 
ery space and classes of extragalactic sources that were poorly 
sampled at X-ray energies in the past (e.g., star-forming and 
passive galaxies, and obscured and low-luminosity AGNs). 

Bauer etal. (2004; hereafter, B04) measured the number 
counts for X-ray point sources detected in the «2 Ms Chandra 
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Deep Field-North (CDF-N) and wl Ms CDF-South (CDF-S), 
the deepest Chandra surveys at the time. Using the available 
multiwavelength data, B04 were able to distinguish between 
different X-ray emitting populations (e.g., active galactic nu- 
clei [AGNs], normal galaxies, and Galactic stars), and mea- 
sure their contributions to the CXRB. B04 found that, gener- 
ally, AGNs make up exclusively the relatively bright number 
counts. However, at the faintest flux levels in the 0.5-2 keV 
bandpass, the most sensitive bandpass studied in the CDF sur- 
veys, B04 showed that normal galaxies start to comprise an 
appreciable fraction of the number counts (~25% of sources 
with 5o.5-2 keV ~ [2-10] X 10"'^ ergs cm'^ s"'). 

Recently, the deepest extragalactic X-ray survey yet con- 
ducted, the CDF-S, has reached a total exposure of «4 Ms 
(Xue etal. 2011; hereafter, Xll). Xll presented point-source 
catalogs and data products for the survey and provided ba- 
sic multiwavelength classifications for the 740 individually 
detected X-ray sources in their main catalog. This work 
has revealed that, at the faintest flux levels, normal galax- 
ies are playing an increasingly important role in the new 
sources detected. For example, at 0.5-2 keV fluxes below 
~5 X 10"'^ ergs cm"^ s"', «50% of the X-ray detected sources 
are classified as likely normal galaxies; below the current de- 
tection limits, it is almost certain that the normal galaxy frac- 
tion continues to increase (e.g., B04; Ranalli etal. 2005). 

In this paper, we present number counts for the new «4 Ms 
CDF-S, focusing on the faint-flux regime (0.5-2 keV fluxes 
^ 10"'^ ergs cm"^ s"'). In §2, we describe a new Bayesian 
method for computing number counts, which properly ac- 
counts for biases and measurement uncertainties that are 
present in the important flux regime near the detection limit. 
In §3, we utilize the available multiwavelength photometric 
data as well as optical spectroscopic and photometric redshift 
catalogs to quantify, with good reliability, the relative contri- 
butions to the number counts from AGNs, normal galaxies, 
and Galactic stars. We further break down the number counts 
to quantify the contributions that subpopulations make to the 
AGN (with redshift, intrinsic AGN column density, and X-ray 
luminosity subpopulations) and normal galaxy (with redshift 
and morphological subpopulations) number counts. We con- 
clude §3 by discussing how each population and subpopula- 
tion contributes to the overall CXRB intensity. In §4, we use 
our subpopulation number-counts estimates to compare with 
predictions from phenomenological models that describe how 
the X-ray emission from accreting SMBHs and normal galax- 
ies are expected to evolve with cosmic time. In §5, we use 
these models to estimate directly number counts to flux levels 
below current detection limits and highlight the prospects of 
deeper X-ray surveys. In §6, we summarize our findings. 

Throughout this paper, we make use of the main point- 
source catalog and data products provided by Xll. The 
Galactic column density for the CDF-S is 8.8 x 10'^ cm~^ 
(Stark etal. 1992). All of the X-ray fluxes and luminosities 
quoted throughout this paper have been corrected for Galac- 
tic absorption. In the X-ray band, we make use of four 
bandpasses: 0.5-2 keV (soft band [SB]), 2-8 keV (hard band 
[HB]), 4-8 keV (ultra-hard band [UHB]), and 0.5-8 keV (fuU 
band [FB]). Values of Hq = 70 km s"' Mpc"', Om = 0.3, and 
Qa = 0.7 are adopted throughout this paper (e.g., Spergel etal. 
2003). 



2. METHODOLOGY FOR COMPUTING NUMBER COUNTS 
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Fig. 1 . — Detection probability as a function of intrinsic 0.5-2 IceV flux for 
point sources located on axis (dashed curve) and f»8' off axis (solid curve). 
These curves were computed considering a counts-to-flux conversion fac- 
tor appropriate for a power-law SED with F = 1.4 (see eqns. 2 and 3), the 
mean photon index of the CXRB. Each curve asymptotes to the defined de- 
tection probability threshold Pihresh = 0.004 (dotted line) at low flux levels. To 
construct the cumulative solid angle as a function of flux limit for the entire 
Ki4 Ms CDF-S image, we computed such curves at each location on the Chan- 
dra image multiplied by the effective soild angle per pixel (RiO.242 arcsec^) 
and added them (see Fig. 2). 

Our primary goal is to evaluate the number counts across 
the entire «465 arcmin^ «4 Ms CDF-S to the faintest possible 
flux levels. In the faint-flux regime, computing number counts 
presents a challenge, since these calculations must properly 
account for (1) the non-negligible spatial variations in sen- 
sitivity across the Chandra image, (2) incompleteness issues 
related to source detection algorithms (see §2.1), and (3) the 
Eddington bias. The B04 investigation of the CDF number 
counts implemented Monte Carlo simulations using an ex- 
trapolated faint-flux number-count model to measure and cor- 
rect for these biases. Although this method provides reason- 
able first-order corrections that account for the biases near the 
flux limits, it does not optimize the input faint-flux model. 
Such an optimization would require several Monte Carlo runs 
with varying faint-flux extrapolations, which is not feasible 
due to the computational reqirements of this procedure. Such 
a limitation will therefore introduce small systematic errors 
in the number-count measurements near the flux hmit if the 
wrong faint-end number-count prior is chosen. 

To mitigate this difficulty, we employ an approach similar 
to the Poissonian-based methods described by Georgakakis 
etal. (2008; hereafter, G08); however, our method has been 
adapted to account for the different source-detection meth- 
ods adopted for the «4 Ms CDF-S catalog by Xll. The 
XI 1 source-detection criteria make use of both wavdetect 
(Freeman etal. 2002) for initial source selection and ACIS 
Extract (AE) to improve photometry and re-evaluate 
source-detection significance (Broos etal. 2010). AE makes 
use of binomial statistics to evaluate source significance, and 
it properly accounts for PSF variations between observations 
and uncertainties in local background measurements (see be- 
low). Our method for computing number counts, described 
in detail below, uses a Bayesian approach with maximum- 
likelihood optimizations to account for the Eddington bias and 
completeness hmitations without requiring a large number of 
time-consuming simulations. 



X-ray Number Counts in the 4 Ms CDF-S 
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2.1. Source Recovery Functions and Flux Probability 
Distributions 

Number-counts computations at the flux Hmits depend sen- 
sitively on the solid angle of the survey over which a source 
of a given flux could be detected. As described in §3 of Xll, 
to be included in the main catalog, an X-ray source must (1) 
be detected by wavdetect at a false-positive probability 
threshold of 10"^ and (2) contain s counts (derived from AE; 
see column 8-16 of Table 3 in Xll) within an aperture rep- 
resenting «90% of the point-source encircled-energy fraction 
(EEFrT that satisfies the following binomial probability cri- 
terion: 

Fix >s) = y] p-'H -p)"-' < Pthresh, (1) 

^-^ xUn-xy. 

where n = s+b^^t and p = 1/(1 +^ext/^si-c)- Here b^^t is the to- 
tal number of background counts extracted from a large region 
outside of the point source (while masking out regions from 
other X-ray detected sources) that was used to obtain an esti- 
mate of the local background count rate. The quantity b^rc is 
the estimated number of background counts within the source 
extraction region, which was measured by rescaling feext to the 
area of the source aperture (typically, /^ext/^src ~ 16; see §4.1 
of Xll for details). In equation 1, Pthresh = 0.004, the value 
adopted by Xll. This choice of Pthresh was empirically se- 
lected both to optimize the number of sources detected and to 
ensure that nearly all detected sources are reliable (see §4. 1 of 
XI 1). This multi-stage procedure for identifying a highly re- 
liable list of source candidates will not result in a complete 
selection of all real sources in the image having binomial 
probabilities < Pthresh, since the initial wavdetect selec- 
tion has more complex source-detection criteria (see Freeman 
et al. 2002 for details) than the simple criterion given in equa- 
tion 1 (see §2.2 and the Appendix for discussion on correcting 
for these issues when computing number counts). 

Using a probabilistic approach, we compute the solid angle 
of the sky within which a source with intrinsic flux S would be 
detected if present. First, using eqn. (1), we define the quan- 
tity L, which is the number of counts required for a detection, 
which satisfies the relation P(x >L) = Pthiesh = 0.004. In order 
to compute L at each location in the CDF-S image, it is nec- 
essary to estimate local values of /jsi-c and /jext that are similar 
to those that would be estimated following the Xll approach 
used for the X-ray detected sources. Unfortunately, the Xll 
procedure for extracting these values at every location in the 
image is computationally prohibitive, since source photome- 
try is performed on an observation-by-observation basis (i.e., 
in up to 54 observations in total) and makes use of polygonal 
apertures that approximate local PSFs (see Xll for further de- 
tails). To overcome this issue, we estimated Z^src by extracting 
background counts from the merged CDF-S background maps 
(see §7.1 of Xll) using circular apertures with sizes that en- 
compass the 90% EEF for a point source. This approach was 
tested by comparing our /jsrc values measured in the regions 
of the Xll sources with those of Xll. We find good agree- 
ment between values and a Icr scatter at the w 18-25% level. 

As noted by XI 1, in special cases where sources liad nearby neighbors, 
the source-detection probability was measured using an aperture smaller than 
the Si90% EEF; however, only ^2-3% of all sources had moderate source 
crowding, where detection probabilities were measured using apertures that 
encompass ^ 50% EEF. We therefore expect that this will have a negligible 
effect on our number-counts measurements. 
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Fig. 2. — Effective solid angle A(S) versus intrinsic flux limit S for the SB 
(dotted curve), HB (short-dashed cur\>e), UHB (dashed-dotted cun'e), and 
FB (solid cun'e). These curves have been computed for a power-law SED 
with r = 1.0, the median photon index for sources detected in both the SB 
and HB. The shaded regions give the effective solid angle curves appropriate 
for the interquartile (f»25-75%) range F = 0.4-1.5. These curves were com- 
puted following the probabilistic methods discussed in §2.1 and account for 
uncertainties in measured flux conversions at the lowest flux levels where the 
number of detected counts are small. The horizontal dotted line represents 
the Si 10 arcmin^ limit, above which we can confidently compute number 
counts (see §2.1 for details). These curves are qualitatively different in na- 
ture from those produced by Xll, which utilize single- valued count-rate to 
flux conversions. 

To estimate appropriate values of /jgxt at each location (pixel) 
on the CDF-S image, we measured the local off-axis angle 9p; 
using the XI 1 point-source catalog, we adopted the maximum 
value of /jext for sources with off-axis angles 9 = 9p± Q'.25. 
Given values of b^rc and b^^t, we numerically solved the rela- 
tion P{x > L) = 0.004 to obtain L at each image location. In 
this manner, we constructed a spatial sensitivity map consist- 
ing of L values. 

In principle, L could be combined with the exposure time 
to estimate a count-rate limit, which can in turn be converted 
to a single flux limit using a count-rate to flux conversion (as 
done in the Xll catalogs); however, this single flux-limit ap- 
proach does not directly incorporate the probabilistic nature of 
source detection that is important in the low-count regime. It 
is therefore more informative to compute the probability that 
a source with intrinsic flux, S, would be detected given that L 
counts are required. Such a source is expected to contribute 
the following number of counts to the source detection cell: 

T = tt,^pCr]S+b^rc (2) 

where fexp, C, and 77 are the effective exposure time, the con- 
version from flux to count-rate, and the encircled-energy frac- 
tion, respectively. Values of fexp are taken directly from the 
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exposure maps from XI 1 (see their §3.1). These exposure 
maps include the effects of vignetting, gaps between CCDs, 
bad-column filtering, bad-pixel filtering, and spatial and time- 
dependent degradation in quantum efficiency due to contami- 
nation on the ACIS optical-blocking filters. Values of C will 
depend on the spectral energy distribution (SED) shape. In 
this work we use power-law SEDs to characterize the ob- 
served spectra; these SEDs can be described using the photon 
index T. For each source, F was derived using the HB-to- 
SB count-rate ratio (corrected for differential vignetting and 
exposure times) as a proxy for spectral slope. The HB-to- 
SB count-rate ratio was calibrated against F using the AE- 
automated XSPEC-fitting procedure for relatively bright X- 
ray sources (with full-band counts greater than 200; this en- 
sures reUable XSPEC-fitting results). This approach takes 
into account the multi-epoch Chandra calibration information 
(see XI 1 do further details). 

Therefore, at each location on the image, the probability 
of source detection as a function of flux S can be computed 
following 



niet 



L+bcx, 

= E 

x=L 



(3) 

where /jext is the extracted background counts that were used 
to calculate the local background in eqn. (1). Figure 1 shows 
/"(let as a function of 0.5-2 ke V flux 5o.5-2 kev for randomly se- 
lected positions at ^ « 0' and 8' offsets from the average aim- 
point of the 4 Ms CDF-S (assuming r/ = 0.9 and a counts-to- 
flux conversion factor appropriate for a power-law SED with 
r = 1 .4). We note that as the flux drops to zero, Pdet asymptot- 
ically approaches our detection threshold Pthi-esh = 0.004. This 
shows that, even when no source is present, there is still a 
finite (yet smaU, i.e., P = 0.004) probability that a positive 
fluctuation may exceed our adopted detection threshold; how- 
ever, as shown by Xll through multiwavelength counterpart 
matching, the initial wavdetect source selection ensures 
that very few false sources are present in the Xll catalog. 

We consider that each Chandra pixel represents a 
smaU local solid angle of size d^l « 0.242 arcsec'^ (i.e., 
«0?'492 X Of'492) with a flux and SED (i.e., V) dependent 
detection probability distribution characterized by fjet- Un- 
der this consideration, the total surveyed solid angle over 
which sources with flux S and photon index F could be de- 
tected when present is therefore A(5,F) = X],/^det,! (here- 
after, the effective solid angle), where the summation is over 
all possible detection cells (i.e., all pixels). In Figure 2, we 
show A(5, F = 1 .0) versus 5 for flie «4 Ms CDF-S in flie four 
bandpasses. The shaded regions show A(5,F) in the range 
of F = 0.4-1.5 (median value F = 1.0), which represents the 
interquartile (i.e., 25-75%) range for the 332 Xll sources 
with estimates of F that were not based on limits. We note 
that generally A(5,r) increases with increasing F. We find 
that the A(5,F) curves asymptote to a value of k,2 arcmin^ 
approaching 5 = 0, suggesting that sources with extremely 
low fluxes (down to zero) could in principle produce fluctu- 
ations exceeding the probability threshold defined in equa- 
tion 1. However, as discussed above, such sources would 
most efficiently be removed using our initial wavdetect 
screening. Therefore, we cannot use information at such flux 
levels to determine number counts reliably. We therefore 
choose to restrict our number counts computations to flux lev- 
els where >10 arcmin^ solid angle is accessible in our survey 
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Fig. 3 . — 0.5-2 keV flux probability distributions for two AGNs in the XI 1 

catalog with no-source detection probabilities (provided by AE) of silO~^^ 
(solid curve) and k().()()3 (dashed curve). The latter source is roughly at the 
boundary of our adopted detection threshold Pihresh = 0.004, and therefore 
illustrates the extent to which the flux conversion for sources in the low-count 
regime with non-neghgible Eddington biases can be affected. The vertical 
dotted lines indicate the most probable single-value flux given in the XI 1 
catalog. 

for the case of F = 1.4 (the mean SED of the X-ray back- 
ground; e.g., Moretti etal. 2009). Our adopted «10 arcmin^ 
solid-angle limit additionally constrains our flux limits. The 
resulting Umits are 5.1 x 10""^, 3.7 x 10"'^ 4.6 x lO"''', 
2.4 x 10"" ergs cm^^ s"^ for the SB, HB, UHB, and FB, re- 
spectively; these are factors of « 1.9-4. 3 times fainter than 
those of B04 and G08, the previously deepest number-counts 
studies. For the full range of X-ray spectral slopes (i.e., F 
values) in the Xll sample, «5-20 arcmin^ of solid angle is 
available for number counts computations at the flux Hmits 
(above the asymptotic regime). We note that the flux limits 
derived here are fainter than those presented in Xll, which 
were «9.1 x 10"'** and «5.5 x 10"'^ ergs cm"^ s"' in the SB 
and HB, respectively (see also the area curves in our Fig. 2 
compared with Fig 23 in XI 1). This is due to the fact that XI 1 
considered only a single count-rate to flux conversion factor 
and did not use the probabilistic approach adopted here. 

To account for the fact that, for each X-ray detected source 
in our main catalog, we are only able to measure reliably the 
total observed counts s (see above), the local background b^xt, 
and that the intrinsic flux S may be subject to large uncer- 
tainty (particularly in the low-count regime), we consider the 
conversion from counts to flux for each source to be prob- 
abilistic. For each X-ray detected source, we computed the 
flux probability distribution foUowing 



(s + beuV- 



ilfcext! \ T + h 



T+be 



dN/dS\ 



model 



(4) 

where the term dN/dS\^^^^^ is a Bayesian prior, based on the 
differential number counts, which accounts for the Eddington 
bias near the sensitivity limit. As noted in B04, the slope of 
the number counts of AGNs, normal galaxies, and Galactic 
stars will differ at the flux limit of the CDF-S. Therefore, 
dN /dS\^^^^^ for a given source will depend on which source 
population it belongs to. Previous studies of X-ray number 
counts (e.g., Rosati etal. 2002; B04; Kim etal. 2007; GOB) 
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Fig. 4. — Observed-frame 0.5-8 keV luminosity versus redshift for AGNs 
(blue triangles) and normal galaxies (red squares). Luminosities were calcu- 
lated using the best redshift estimates and 0.5-8 keV fluxes provided by Xll 
and our adopted cosmology. The dashed curve indicates the observed-frame 
luminosity corresponding to the flux limit at the center of the Ri4 Ms CDF-S. 
The normal-galaxy popluation broadly covers the redshift range z ~ 0.03-2.6 
and AGNs cover the redshift range of z ~ 0.1-8. These redshifts are based 
on both spectroscopic and photometric redshifts, and all sources at z > 4.76 
have only photometric redshift estimates. We adopt the most probable pho- 
tometric redshifts; however, some of the z > 4.76 photometric redshifts are 
consistent with being at lower redshifts (see Luo et al. 2010 for details). 

have shown that power laws provide good fits to the overall 
shapes of the logA^-logS'. To first order, we use priors based 
on the following power-law parameterizations: 



dN 

Is 



AGN 



=/:°'''(5/5,ef)- 



'dS 

J »7 Star 

dS 



(5</r^) 



(5) 



where dN/dS{AGN), dN/dS(gci\), and dN/dS(stSii) are dif- 
ferential number-count parameterizations to be applied to 
AGNs, normal galaxies, and Galactic stars, respectively, 
f^'^^ is the flux related to the break in the double power- 
law used to describe the AGN number counts, and 5,ef = 
10"'"^ ergs cm"^ s"'. As we show in §3.1 below, we char- 
acterize each X-ray source using the X-ray and multiwave- 
length data and provide best estimates of the parameters in 
equation 5 (i.e., K, f3, and fb values) for the AGNs, nor- 
mal galaxies, and Galactic stars. For each source, we com- 
puted Pj down to the flux limits defined above and normalize 
equation 4 using J Pj dS = 1. In Figure 3, we show exam- 
ples of Pr as a function of 0.5-2 keV flux for two AGNs in 
the main Chandra catalog having AE probability P w 10""^^ 
(solid curve) and P ss 0.003 (dashed curve; near the detection 
threshold). For this computation, we utilized values of (3^'^^ = 



1.49,/?^™ = 2.49, and/b'^™' 



: 5.6 X 10 ergs cm ^ s ' (the 



normalization is arbitrary); as we will show in §3.1, these val- 
ues represent the best-fit parameterizations for AGN SB num- 
ber counts. 

2.2. Cumulative Number Counts Computation 



Using the above information, we computed cumulative 
number counts for the w4 Ms CDF-S survey using the fol- 
lowing integral: 



N(>S) = 



J2PT.i(S)C'i/A(S,Td 



dS, 



(6) 



where the summation is over all X-ray detected sources in a 
particular category of sources (e.g., all sources, AGNs, nor- 
mal galaxies, etc.) and F,- represents the effective photon in- 
dex for source /. For the 332 X-ray sources detected in both 
the SB and HB, we utilized the photon indices from XI 1 ; oth- 
erwise, we used F = 1 .4. Since Pj j is dependent on our input 
Bayesian prior (which depends on source type), the number 
counts will depend mildly on the choice of our model. The 
term C' is a count-rate and off-axis angle dependent complete- 
ness term, which corrects for the fact that our X-ray point- 
source catalogs were constructed by (1) running wavdetect 
to form an initial list of candidate sources and (2) assessing 
the probability of detection using AE. If we chose to use a 
catalog based strictly on the AE source-detection probability 
assessment for every location on the image, then C' = 1 for all 
sources; however, such a catalog would produce large num- 
bers of false sources, leading to significant errors in the num- 
ber counts at faint fluxes, and is computationally impractical. 
In practice, the completeness corrections are small perturba- 
tions (median C' « 1.00 for all bands, with «95% of sources 
having C' ^ 1.1-2 A) and affect number counts near the flux 
limits. C' and its count-rate and off-axis angle dependencies 
were measured using marx simulations, and a full description 
of this procedure can be found in the appendix. 

Following the same approach implemented in §5 of G08, 
we utilized maximum-likelihood techniques to optimize the 
Bayesian model parameters. As discussed above, we made 
use of power-law differential number-counts models for our 
Bayesian priors (see equation 5), which are dependent on 
source type. For a given source, we first characterized its 
type (i.e., AGN, normal galaxy, or Galactic star; see §3.1 be- 
low), chose an appropriate model in equation 5, and computed 
the probability of that source being present in the «4 Ms 
CDF-S catalog following 



J PT(S,Ni)dS 



IdN/dSl^,^,A(S,r,)dS' 



(7) 



where A^, represents the total number of counts measured 
within the extraction cell for each detected source. Consid- 
ering our model, the total likelihood of obtaining the «4 Ms 
source catalog and its source-count distribution can be com- 
puted as YiiPi- We maximized the total likelihood for the 
model input parameters for AGNs, normal galaxies, and 
Galactic stars separately to find best-fit values for each of the 
parameters in equation 5. 

3. RESULTS 

3.L Number Counts by Source Type and Total Counts 

The extensive multiwavelength data in the «4 Ms CDF-S 
region allow for the robust characterization of the X-ray de- 
tected sources. As described in Xll, 716 of the 740 X-ray 
detected sources have multiwavelength counterparts, and 673 
have either secure spectroscopic or reliable photometric red- 
shifts. Following a similar scheme to that provided by Xll, 
we have classified the X-ray sources in our sample as active 
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Fig. 5. — Top panels of a-d: Cumulative number counts for the SB (a), HB (b), UHB (c), and FB (rf) broken down into AGNs (open blue triangles), normal 
galaxies {open red squares), and Galactic stars (open green stars). Total number counts have been shown as filled circles. In each plot, the best-fit dN /dS 
parameteiizations based on equation 5 have been shown as blue long-dashed, red dotted, and green short-dashed curves for AGNs, normal galaxies, and Galactic 
stars, respectively. The total number counts model, based on summing the three contributing components, has been shown as a black curve. Bottom Panels of 
a-d: Fractional contributions from AGNs, normal galaxies, and Galactic stars to the total number counts. For the majority of the flux ranges, AGNs dominate the 
number counts; however, nornial galaxies provide significant contributions near the flux limits of the SB and FB. [A machine-readable table of the cumulative 
number counts data for the ^4 Ms CDF-S is provided in the electronic edition.] 
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Fig. 6. — Differential number counts (dN/dS) versus flux for the SB (a), HB (h), UHB (c), and FB (d), broken down into AGNs (open blue triangles), normal 
galaxies {open red squares), and Galactic stars {open green stars; for the SB, HB, and FB). The differential number counts have been estimated using flux bins 
of faO.3 dex. In each plot, the best-fit dN jdS parameterizations based on equation 5 have been shown as blue long-dashed, red dotted, and green short-dashed 
curves for AGNs, normal galaxies, and Galactic stars, respectively. The total number counts model, based on summing the three contributing components, has 
been shown as a black curve. Bottom Panels of a-d: Fractional contributions from AGNs, normal galaxies, and Galactic stars to the differential number counts. 
In the SB and FB, we find that the differential number counts of normal galaxies appear to surpass those of AGNs for fluxes below ^(1-2) X 10"" ergs cm"- s"' 
and fs(3-7) X 10"" ergs cm"- s"', respectively. 
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TABLE 1 

MAXIMUM-LlKIiLlllOODBliS r-Fl r MODliL Paraml i lrs 







AGNs 






Normal Galaxies 


Galactic Stars 


Bandpass 
(1) 


^AGN 

(2) 


^AGN 

(3) 


^AGN 

(4) 


rAGN 
■'break 

(5) 


^gal 

(6) 


^gal 

(7) 


t-slar 

(8) 


^star 

(9) 


0.5-2 keV 

2-8 keV 

4-8 keV 

0.5-8 keV 


169.56 ± 8.69 
573.13 ±27.49 
1463.27 ± 63.61 
562.20 ± 22.96 


1.49 ±0.03 
1.32±0.04 
1.08 ±0.10 

1 ^440.04 


2.48 ±0.27 

2 55+0" 
2.23 ±0.09 
2.35 ±0.15 


6.4±1.0 
8 1+1-5 


1.53 ±0.10 
1.10 ±0.22 
0.27 ±0.10 
2.82 ± 0.26 


r, O^+U.08 

2.29 ±0.25 

2 4^+0.38 


3.97 ± 0.37 
0.64 ± 0.22 

4.07 ± 0.51 


1 41+0.14 
^•^^-0.13 

1.79 ±0.50 

1 55+0- 19 
^••'■'-n.is 



NOTE. — Best-fit values and la errors for our number-counts priors. Column (1) lists the bandpass. Columns (2)-(5) provide the double power-law differential 
number-counts parameterization for AGNs, including the normalization K^^^ (Col. 2; in units of lO''' deg"^ [ergs cm"^ s"']"'), faint-end slope (Col. 3), 
bright-end slope ^2'^^ C^^l. 4), and break-flux Z,^^ (Col. 5; in units of 10"!^ ergs cm~^ s"'). Columns (6) and (7) provide the normal-galaxy single-power-law 

normalization and slope respectively. Columns (8) and (9) provide the Galactic star single-power-law normalization Kf^ and slope respectively. 
These values were computed using the maximum-likelihood methods described in §2.2. 



galactic nuclei (AGNs), normal star-forming galaxies (here- 
after, galaxies), and Galactic stars. In this process, we make 
use of six AGN selection criteria: 

1. X-ray Luminosity: A source with an observed-frame 
SB, HB, UHB, or FB luminosity (i.e., observed- 
frame Lx = Andlfx; where is the luminosity dis- 
tance) ^ 3 X 10''^ ergs s"' is classified as a luminous 
AGN. In the nearby universe, the most powerful star- 
forming galaxies (luminous infrared galaxies) all have 
io.5-8 keV ergs s"' (e.g., Iwasawa etal. 2009, 
2011; Lehmeret al. 2010; Pereira-Santaella etal. 2011). 



2. X-ray Spectral Shape: A source with an effective pho- 
ton index of F ^ 1 .0 is most hkely to be dominated 
by a single powerful X-ray source that is significantly 
obscured. We consider such sources to be obscured 
AGNs. 

3. X-ray to Optical Flux Ratio: A source with an X-ray 
to optical flux ratio of log(/x//x) > 1 (where fx = 
/0.5-2 kev, fi'S kev, /4-8 kev, or /o.5_8 kev) IS Considered 
to have X-ray emission significantly elevated compared 
with normal-galaxy processes (e.g., hot gas, low-mass 
and high-mass X-ray binary emission) as traced by stel- 
lar emission (i.e., /r). Such sources were classified as 
AGNs. 



4. X-ray to Radio Luminosity Ratio: A source with sub- 
stantial excess (i.e., a factor of ^ 5) X-ray emission 
over the level expected from pure star formation (i.e., 
observed-frame SB, HB, FB, or UHB luminosity ^ 5 x 
[8.9 X lO'^Li.4 GHz]; Bauer etal. 2002) is considered to 
be X-ray overluminous compared with the X-ray/SFR 
correlation and its X-ray emission is therefore likely to 
be dominated by an AGN. This criterion is similar in 
nature to that of criterion 3, but differs in that it allows 
for a more sensitive classification of a subset of actively 
star-forming galaxies (i.e., detected at 1.4 GHz) using a 
more reliable tracer of the intrinsic galactic SFR that is 
insensitive to dust obscuration. 

5. X-ray Variability: Sources found to be variable in the 
X-ray band (with probability ;^ 95%) on timescales of 
months to years that also have 0.5-8 keV luminosities 
greater than 10*^ ergs s~^ are unlikely to be produced 
by normal-galaxy processes (see Young etal. 2012 for 
details). These sources were classified as AGNs. 

6. Optical Spectroscopy: Finally, sources with optical 
spectroscopic AGN features such as broad and/or high- 
excitation emission lines are classified as AGNs. These 
sources were identified using the spectroscopic catalogs 
of Szokoly etal. (2004), Mignoh etal. (2005), and Sil- 
verman etal. (2010) (see XI 1 for details). 



TABLE 2 
Number Counts Statistics 



Number 



Total 
(deg-2) 



CXRB Fraction 
(%) 



Class 


Subclass 


SB 


HE 


UHB 


FB 


SB 


HB 


UHB 


FB 


SB 


HB 


UHB 


FB 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(1) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


-T Dngni-v_-oiTecuoii 


lOLdl 


















13.1 m J 


9.1 A A- 1 A 


9.9. A A~ \ 9 
oo.'+ IL U.o 


CI f. A- 9 Q 

o 1 .o IC o.y 


AU 


Total 


650 


403 


260 


634 


27832 ± 1803 


10495 ± 871 


8387 ± 787 


22579 ± 1506 


53.1 ±3.4 


66.3 ± 5.8 


64.6 ± 6.7 


62.7 ± 4.2 


AGNs 


Total 


474 


387 


256 


520 


14925 ± 1228 


9310 ± 776 


8053 ± 757 


12802 ± 943 


49.0 ± 4.0 


65.5 ± 5.8 


64.1 ± 6.6 


59.6 ± 4.4 




uiikiiowii z 
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AO 
'fZ 


ZDWJ HI J^J 


lujo m jZD 


f,00 A- OAQ 

ozz ^ zf y 


lOflA -t- A19 


1.1 m w.z 


1 Q -1- n ^; 
1.7 m u.D 


1 Q -1- n R 
i.y m u.o 


A 9. A- (\ A 
l.o IC U.f 




z = 0.0-1.5 


209 


182 


123 


238 


5630 ±721 


3929 ± 478 


3677 ± 537 


5077 ± 523 


3 1.6 ±4.0 


40.0 ± 5.2 


37.9 ±6.1 


37.3 ±3.8 




z= 1.5-3.0 


165 


133 


89 


181 


5338 ± 763 


3111 ±428 


2780 ±413 


4738 ±612 


12.9 ± 1.8 


17.7 ± 2.6 


17.3 ±2.8 


15.9 ±2.1 




z> 3.0 


57 


51 


35 


59 


1353 ± 328 


1213 ±289 


975 ± 228 


1082 ± 266 


3.4 ±0.8 


5.9 ± 1.5 


7.1 ± 1.8 


4.7 ± 1.1 




logWH < 22.0 


102 


86 


51 


108 


993 ± 105 


1559 ± 262 


828 ± 159 


1572 ± 40 


30.3 ± 3.2 


21.5 ± 3.9 


16.1 ± 3.4 


24.7 ± 0.6 
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Zf 7 


Do / I III OlO 


III J JO 


A\ ^7 -1- 

H-l J / III JOJ 


OJ't-o III 01 


\A A A- \ 1 
I't.'t ^ 1. / 


9Q n -1- ^ 7 


9S 1 A- A'X 
Zo.l III *+.J 


OA "^J^ A 'X 




logiVn > 23.0 


105 


103 


81 


121 


4458 ± 728 


2105 ± 339 


2447 lb 382 


2777 lb 54 


3.1 ±0.5 


13.0 lb 2.2 


18.0 ± 3.0 


8.9 4 0.2 








zo 


1 9 

IZ 




ZZoj It J't4- 


007 III Z'fO 




9111 -1- AO'X 
Z 1 1 1 ^ f Z J 


\ f\ A- (\ A 
l.D III U.'t 


1 7 -1- n s 
1. / m U.J 


Z.U III U.o 


1 Q & A A 

i.y at u.f 




logZ-x = 42.0-43.0 


124 


86 


40 


142 


5086 ± 704 


3247 ±514 


1592 ± 359 


5348 ± 677 


3.4 lb 0.5 


6.3 ± 1.0 


6.1 ± 1.5 


5.3 §.0.7 




log Lx =43.0^4.0 


167 


159 


113 


186 


3818 ±626 


3067 ± 388 


3611 ±485 


2534 ± 267 


12.4 lb 2.0 


21.6 ± 2.9 


24.4 ± 3.6 


17.99i 1.9 




Inn T -.r /I /I A 


yo 


OS 


S9 
oZ 


1 ni 


1 1 a/1 _l_ 1 /;7 

1 1 ^ lO / 




i^ZJ ^ Zj^ 


QC\A -1- 1 CO 
VU'-f ^ lUZ 


A A- A ^ 


'^A A -1- S 9 
j'4-.U ^ J.Z 


ZV.O ^ J. 4- 


iC^ 7 

c: 


Nonnal Galaxies 


Total 


166 


14 


4 


106 


12704 ± 1316 


1142 ± 394 


334 ±216 


9638 ± 1172 


3.0 ± 0.3 


0.8 ± 0.3 


0.5 ± 0.3 


3 

2.6 £ 0.3 




z < 0.6 


82 


11 


3 


63 


4924 ± 792 


975 ± 371 


165 ± 155 


5049 ± 825 


1.6 ± 0.2 


0.7 ± 0.2 


0.3 ± 0.3 


5' 

1.6^0.3 




0.6 


84 


3 


1 


43 


7770 ± 1052 


167 ± 133 




4573 ± 831 


1.4 ± 0.2 


0.1 ±0.1 




1.0 1^0.2 




Late-Type (star-forming) 


123 


10 


1 


74 


10174 ± 1191 


1051 ± 390 




7222 ± 1014 


2.2 ± 0.3 


0.6 ± 0.2 




1.8^0.3 




Early-Type (passive) 


43 


4 


3 


32 


2530 ± 561 


91 ±53 


294 ± 210 


2416 ± 588 


0.8 ± 0.2 


0.2 ±0.1 


0.4 ± 0.3 


0.8 g 0.2 


Stars 


.... Total 


10 


2 




8 


202 ± 109 


34 ±32 




139 ± 88 


1.2 ± 0.6 


0.1 ±0.1 




0.5 0.3 



"CXRB fractions have been computed using the number-counts measurements presented in this paper below 5 X 10"'', 1.4 X 10"''', 8.4 X 10~", and 2 X 10~" ergs cm~^ s~' for the SB, HB, UHB, and FB, 
respectively, plus bright-end corrections from the Kim etal. (2007) number-counts relations (see § 3.4 for details). 

NOTE. — Number-counts statistics for the source classifications (Column 1) and subclassifications (Column 2) discussed in this paper Columns (3)-(6) provide the number of sources in each subclass used in 



calculating number counts for the four bandpasses. Columns (7)-(10) provide the number counts for each subcategory at the limiting fluxes for the SB (5.1 X 10 ergs cm " : 



'), HB (3.7 X Ur" ergscm-2 . 



'), 



UHB (4.6 X 10" 



s ), and FB (2.4 X 10 ergs cm s ), respectively. Columns (11)-(14) provide the integrated intensity (expressed as the percentage of the CXRB intensity) for the SB, HB, UHB, 



andFB, assuming CXRB intensity values of (8. 15 ±0.58) X lO''^, (1.73 ±0.23) x 10"", (1.04 ±0.14) x 10"", and (2.54 ±0.24) x 10"" ergs cm'^ s"' deg'^, respectively (see § 3.4 for details). 
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Fig. 8. — Euclidean-normalized, cumulative number counts (A'[> S] X [S/5,cf]'^, where 5i-ef = 10"''' ergs cm"^ s"') and Icr errors for the Si4 Ms CDF-S 
{black cun'es with gray error bars) for the SB (a), HB (h), UHB (c), and FB {d). The Miyaji & Griffiths (2002) predicted number-counts boundaries (from 
fal Ms fluctuation analyses) are shown as the dotted regions in panels a and h. Previous Ms CDF-N plus Ril Ms CDF-S number counts estimates from B04 
and si2 Ms CDF-S number counts from LOS are indicated as blue dotted and red short-dashed curves, respectively. The G08 number counts, which include the 
Ms CDF-N and wl Ms CDF-S surveys, are indicated as dashed triple-dotted curves. Similarly, the ChaMP survey number counts from Kim etal. (2007), 
which include the Si2 Ms CDF-N and sil Ms CDF-S surveys, have been shown as green long-dashed curves. Number-counts measurements from other XMM- 
Newton and Chandra surveys have been indicated, including the XMM-Newton serendipitous (Mateos etal. 2008; magenta triangles), Lockman hole (Brunner 
et al. 2008; brown squares), and COSMOS (Cappelluti et al. 2009; orange circles) surveys, as well as the Chandra COSMOS (Elvis et al. 2009; see also Puccetti 
et al. 2009; gold stars). In the case of the UHB, we highlight the f» 1 Ms CDF-S number counts from Rosati et al. (2002; cyan dot-dashed curiae). We note that all 
measurements shown (i.e., both curves and data points) correspond to measurements, and models are not plotted here. We find good overall agreement between 
our number counts measurements and those found in previous studies. [A machine-readable table of the cumulative number counts data for the Ri4 Ms CDF-S is 
provided in the electronic edition.] 



X-ray Number Counts in the 4 Ms CDF-S 



11 



A source is initially classified as an AGN if at least one 
of the above six criteria are satisfied. However, we found 
that five of the sources that were classified as AGNs via the 
above criteria had 1-2 keV and 0.5-1 keV band ratios that 
were consistent with being dominated by hot interstellar gas 
(see Danielson etal. 2012). The initial AGN classification 
for these five sources was based on having high X-ray to 
optical flux ratios (criterion 3); however, their X-ray lumi- 
nosities compared with their B-band luminosities appear to 
be consistent with hot-gas dominated galaxies (see Daniel- 
son etal. 2012 for details). We therefore classified these 
five objects as normal galaxies. In addition to classifying 
sources directly as AGNs and normal galaxies, we searched 
for direct multiwavelength indicators of Galactic stars. Likely 
Galactic stars were identified using (1) optical spectroscopy 
from Szokoly etal. (2004), MignoU etal. (2005), and Silver- 
man etal. (2010); (2) HST stellarity indices >0.7 (from the 
Caldwell et al. 2008 GEMS catalogs); and (3) best-fit stel- 
lar templates from MUS YC photometric-redshift fits (Carda- 
mone etal. 2010). All Galactic star candidates were visually 
screened and obvious normal galaxies were noted. The re- 
maining stellar candidates were then classified as "Galactic 
stars." Remaining sources that were not classified as either 
AGNs or stars were classified as normal galaxies. Based on 
these criteria, we therefore estimate that of the 740 X-ray de- 
tected sources, 561 are AGNs, 169 are normal galaxies, and 
10 are Galactic stars. 

We compared our normal-galaxy classifications with the 3 1 
radio and X-ray selected normal galaxies classified by Vat- 
takunnel etal. (2012) and found that 25 sources overlap. We 
find that 4 of the Vattakunnel etal. (2012) sources that we 
classify as AGNs had high X-ray to radio luminosity ratios 
(criterion 4), a criterion that was not implemented by Vat- 
takunnel etal. (2012). 

In Figure 4, we plot the observed-frame 0.5-8 keV lumi- 
nosity versus redshift for the AGNs and normal galaxies with 
redshift measurements. Our AGN and normal galaxy samples 
span redshift ranges of z w 0.1-8 and z ~ 0.03-2.6, respec- 
tively. These redshifts are from both photometric and spec- 
troscopic redshift estimates, and all redshift estimates above 
z = 4.76 are based on photometric redshifts. We have chosen 
to use the most probable photometric redshift for each source; 
however, some of the z > 4.76 sources have photometric red- 
shifts consistent with being at lower redshifts (see Luo et al. 
2010 for details). 

Using the classifications adopted above and our methods 
for computing and optimizing number counts described in §2, 
we computed the cumulative number counts for AGNs, nor- 
mal galaxies, and Galactic stars, and determined the best-fit 
parameters related to their priors (i.e., those characterized in 
equation 5). In Figure 5, we present the breakdown of the 
cumulative number counts from AGNs, normal galaxies, and 
Galactic stars for the four bandpasses (discrete symbols). To- 
tal number counts for the survey are simply the sum of the 
number counts from each of these populations. In Figure 6, 
we show the differential number counts (dN/dS) for the four 
bandpasses. The best-fit dN/dS model parameters, used in 
computing our number counts, have been tabulated in Ta- 
ble 1, and the cumulative and differential number counts cor- 
responding to these models have been shown in Figures 5 and 
6, respectively (continuous curves). In Figure 7, we show 
sample probability distributions for the best-fit dN/dS model 
parameters appropriate for AGNs detected in the SB. 



The bottom panels of Figures 5 and 6 show the fractional 
contribution that the respective cumulative and differential 
number counts of each source type make to the total counts. 
We find that AGNs still dominate the cumulative number 
counts at all flux levels in all four bands; however, at SB and 
FB fluxes below ^ 10"'^ ergs cm~^ s~', normal galaxies un- 
dergo a rapid increase in their numbers. At the faintest SB 
flux limit («5.1 x 10"''^ ergs cm"^ s"'), AGNs and normal 
galaxies reach source densities of w 14,900 and 12,700 deg"^, 
respectively. These sky densities are factors of «2.2 and «4.5 
times larger than those reported at the w2 Ms Chandra depth 
(see B04) for AGNs and normal galaxies, respectively. We 
find that our SB normal galaxy number counts lie between 
the "optimistic" and "pessimistic" estimates from B04 and 
are in good agreement with previous work by Hornschemeier 
etal. (2003). Due to the small solid angle of our survey, the 
bright normal galaxy number counts data drop off at SB fluxes 
4x 10"'^ ergs cm"^ s"^; however, the extension of our best- 
fit power-law fit (red dotted curve in Figure 5fl) agrees well 
with bright-end normal galaxy number counts from the XMM- 
Newton serendipitous plus needles in the haystack surveys 
(Georgakakis et al. 2006) and the Chandra Multiwavelength 
Project (Kim etal. 2006). Normal galaxies comprise 46 ± 
5% and 43 ± 5% of the total cumulative number counts in the 
SB and FB, respectively, a result consistent with that found by 
XI 1, who showed that w40-50% of the SB-detected sources 
within «4' of the w4 Ms CDF-S center are normal galax- 
ies (see their Fig. 13c). This signifies more than a doubling 
in the fractional contribution that the normal-galaxy cumula- 
tive number counts make to the total number counts at the 
SB flux Umit over that found in the ss2 Ms CDFs (see B04). 
In terms of the differential number counts, it appears that the 
contribution from normal galaxies exceeds that of AGNs at 
SB and FB fluxes below w(l-2) xlO~'^ ergs cm"^ s"^ and 
w(3-7) X 10"'^ ergs cm"^ s"', respectively (see Fig. 6a). We 
estimate that near the SB and FB flux limits, normal galaxies 
make up w45-70% of the differential number counts. The 
rapid increase in normal-galaxy counts signifies that deep 
Chandra surveys are on the verge of being normal-galaxy 
dominated at the faintest fluxes (see §5 for a discussion). In 
Table 2, we summarize properties of the number counts. 

In Figure 8, we show our derived total number-count 
distributions normalized by a Euclidean slope of 5''^ (i.e., 
N[> S] X [5/5ref]' ■^ where 5,-ef = 10-'"^ ergs cm'^ s"'). The 
Euclidean representation allows for comparisons with previ- 
ous investigations at both the faint and bright ends of the 
distributions. For comparison, we have plotted number- 
count distributions from previous surveys with XMM-Newton 
(Brunner etal. 2008; Mateos etal. 2008; Cappefluti etal. 

2009) and Chandra (Bauer et al. 2004; Kim et al. 2007; G08; 
Luo etal. 2008; Elvis etal. 2009; Puccetti etal. 2009). Some 
small differences in flux measurements (at the ^ 10% level) 
are expected due to both changes in calibration over the years 
for Chandra\^ as well as cross-calibration uncertainties be- 
tween Chandra and XMM-Newton (see, e.g., Nevalainen et al. 

2010) . In general, we find good agreement between our num- 

For example, as noted in XI 1, in January, 2009 the Chan- 
dra ACIS-I ancillary response file (ARF) was updated resulting in 
a flat Si9% reduction in elfective area and a «:iO-8% reduction be- 
tween 2-5 keV (see http://cxc.harvard.edU/ciao/why/caldb4.l.l_hrma.html I. 
Further uncertainties, at the silO% level, are expected due to spa- 
tial variations of contamination on the optical-blocking filter (see 
[http ://web.mit.edu/iachec/meetings/201 1/Presentations/Marshall.pdf I. 
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Fig . 9 . — Top Panels: Cumulative AGN number counts {filled black circles) in tlie SB (left) and HB (right) broken down by contributions from AGNs in redshift 
ranges of z < 1.5 (open red circles), 7, = 1.5-3 (open green triangles), z > 3 (open blue squares), and no measured redshift (filled gray circles). AGN number-count 
predictions from the AGN population-synthesis models by GOV have been provided both with (dotted curx'es) and without (solid curx'es) an exponential decline 
in the luminosity function at z > 2.7. G07 AGN number-counts models for all redshifts, z < 1.5, z = 1.5-3, and z > 3 have been shown as black, red, green, and 
blue curves, respectively. These models are discussed in detail in § 4. 1 . Middle Panels: Fractional contribution that AGNs in each redshift range provide to the 
overall AGN number counts. The solid curves are the expectations from the GOV model without a declining XLF at z > 2.7. Bottom Panels: The ratio of the AGN 
number-counts data presented in the top panel divided by the G07 model (without a declining XLF at z > 2.7). [A machine-readable table of the redshift-divided 
AGN cumulative number counts dataft>r all four bandpasses is provided in the electronic edition.} 



ber counts and those of previous studies with some minor 
differences. At faint SB fluxes, we find that our number- 
counts measurements are in good agreement with those of 
G08 and B04, which are both based on the w2 Ms CDF-N 
and «1 Ms CDF-S; however, our estimates are w30^0% 
lower than those near the Hmit of the w2 Ms CDF-S (LOS). 
It is likely that this discrepancy is due to an incompleteness 
overcorrection of the faint-end counts in LOS. In the HB, we 
find that our number counts are in good agreement with those 
of B04, but are somewhat lower than those of GOS and higher 
than those of LOS (at the «10-20% level in both cases). In the 
UHB, we find clear evidence of a break in the number counts 
at /4_g kev ~ 10"'^ ergs cm"^ s"' (see Figs. 5-6). Previous 
investigations of the UHB in the CDF-S revealed no obvious 
break down to «7 x 10"'* ergs cm"^ s"' (Rosati etal. 2002), 
which is very close to where we observe the break. We note 
that the GOS investigation found a break at a much brighter 
flux of w6 X 10"'^ ergs cm"^ s"' . In the CDF-S, there are only 
?sl7 UHB sources brighter than this flux, which would make 
it difficult to identify such a break if present. In contrast, the 
GOS analysis made use of wide-area surveys like the EGS, 
ELAIS-Nl, and XBOOTES, and therefore have better source 
statistics at bright UHB fluxes. In all four bands, the «4 Ms 
CDF-S number counts measurements have large uncertainties 
at bright fluxes due to the relatively small solid angle of the 
survey; however, we find good general agreement between 
our measurements and those of previous investigations. 

In §§ 3.2 and 3.3 below, we further divide our number- 
counts measurements into contributions from AGN and 
normal-galaxy subpopulations. For ease of presentation, we 



discuss only results from the SB and HB; however, full re- 
sults from all four bandpasses are provided in the tables for 
reference. 

3.2. AGN Number Counts by Redshift, Intrinsic Column 
Density, and X-ray Luminosity 

As discussed in §3.1, AGNs are the majority of the X-ray 
detected sources and dominate the cumulative number counts 
over all fluxes. The ultradeep Chandra data and redshift in- 
formation allow for first-order estimates of the intrinsic ab- 
sorption column densities, A^h, and 0.5-S keV luminosities, 
Lx, of the detected AGNs. Values of A^h and Lx were esti- 
mated following the procedure outlined in §3.4 of Xue etal. 
(2010). The observed X-ray spectrum was modeled using 
xspec (Amaud 1996) and an assumed absorbed power-law 
model (zpow x wabs x zwabs). The Galactic column den- 
sity was set to S.8 x lO'^ cm"^ (see §1), the intrinsic pho- 
ton index was fixed at Tint = l.S (e.g., Tozzi etal. 2006), and 
the best redshift estimate provided in the XI 1 catalog was 
adopted. Using this model, we determined the A^h value for 
each AGN that reproduced the observed ratio of count rates 
between the HB and SB. For sources detected only in the 
FB an observed spectral index of rob,s = 1.4 was adopted. 
This procedure will be effective for determining A^h and Lx 
values for AGNs in the Compton-thin/relatively unobscured 
regime (i.e., with A^h ^ lO^-'-lO^'* cm"^), and is less reliable 
for a minority of sources that are expected to be Compton 
thick. In the Compton-thick regime, AGN spectra are sig- 
nificantly complex due to the combination of absorption, re- 
flection, scattering, and line emission features that depend on 
metalHcity, ionization state, and geometry (see, e.g.. Murphy 
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Fig. 10. — Cumulative SB AGN number counts for sources at redshifts 
z The Ri4 Ms CDF-S estimates and Itr error bars have been shown 
as open blue squares with error bars. If we add to the z AGN number 
counts the number counts for X-ray detected sources without known red- 
shifts, we obtain the open gray diamonds with error bars. Due to their lack of 
multiwavelength counterparts, in regions where extensive multiwavelength 
coverage is available, the sources with unknown redshifts are good candi- 
dates for being at z ^ 3. Recent constraints from Fiore etal. (2011) for the 
Si4 Ms CDF-S (filled black circles), Civano etal. (2011) for Chandra COS- 
MOS (filled green triangles), and Brusa etal. (2009) fov XMM-Newton COS- 
MOS (filled red squares) have been shown for comparison. The G07 models 
for z?i^ number counts have been shown for the cases where a high-redshift 
(z > 2.7) space-density decline is (blue dotted curve) and is not (blue solid 
curve) included. The recent number-counts predictions, based on the Aird 
etal. (2010) XLF measurement in the 2-10 keV band, ai'e shown as a brown 
long-dashed curve. 

& Yaqoob 2009). Since most of our sources are in the low- 
count regime (i.e., ^ 100 counts), such meaningful modeling 
is beyond our capabilities. Using our power-law model, es- 
timated A^H, and intrinsic flux fx.m estimates, we computed 
Lx = 47rii^/x,int(l +zf '""~^ for each AGN that had a measured 
redshift (spectroscopic or photometric). 

In Figure 9, we display the breakdown of the AGN number 
counts in bins of redshift for the SB {left) and HB (right), 
with the middle panels showing the fractional contribution 
that sources in each redshift range make to the total AGN 
number counts. At the brightest fluxes, low-redshift (z < 1 .5) 
AGNs dominate; however, going to fainter fluxes, we observe 
increasing contributions from z = 1.5-3 AGNs. At the flux 
limit of our survey, AGNs at z = 1 .5-3 make up comparable 
contributions to the number counts as the z < 1.5 AGNs (i.e., 
«40%). For AGNs at z > 3, we find a roughly steady frac- 
tional contribution to the number counts (ss5-15%) across 
the majority of the SB and HB flux ranges where they are 
detected. In Figure 10, we highlight the SB z > 3 AGN num- 
ber counts and compare them with recent studies from the 
literature. We find that our number counts measurements for 
z > 3 AGNs are typically a factor of w 1 .5-2 higher than those 
found in the COSMOS survey fields (i.e., Brusa etal. 2009; 
Civano etal. 2011) and in good agreement with those from 
the early w4 Ms CDF-S study by Fiore et al. (201 1). The dif- 
ference between our values and those of the COSMOS studies 
is likely due to differences in z > 3 AGN selection techniques, 
variations in multiwavelength depth and photometric redshift 



completenesses, and field-to-field variance. 

In Figures 11 and 12, we show the breakdown of AGN 
number counts in bins of intrinsic absorption column den- 
sity, A^H, and 0.5-8 keV luminosity, Lx, for the SB and HB. 
In both bandpasses, we find that unobscured X-ray luminous 
sources dominate the bright-end number counts. Going to 
fainter fluxes, we find that more obscured and less X-ray lumi- 
nous sources make larger contributions. At the SB flux limit, 
AGNs with A^H ^ 10^^ cm"^ and Lx ^ lO'*-' ergs s"' make up 
the majority contributions to the number counts. At the HB 
flux limit, we similarly find that AGNs with A^h 10^~ cm"^ 
dominate the number counts; however, in terms of luminos- 
ity, AGNs with Lx ^ lO'*-' ergs s"' and ^ lO"*^ ergs s"' make 
comparable contributions, with the latter population increas- 
ing faster. 

3.3. Normal-Galaxy Counts by Redshift and Morphology 

As noted in § 3.1, normal galaxies comprise a significant 
fraction (?s32^6%) of the source counts at the SB and FB 
flux limits. X-ray emission from normal galaxies is expected 
to be produced by a variety of populations that differ for dif- 
ferent galaxy populations (see, e.g., Fabbiano 1989, 2006 
for a review). For late-type star-forming galaxies, popula- 
tions associated with a combination of young and old stellar 
populations are expected to provide dominant contributions 
(i.e., high-mass and low-mass X-ray binaries, supemovae and 
their remnants, hot gas from starburst flows, and young stars; 
Colbert etal. 2004; Iwasawa etal. 2009, 2011; Lehmer etal. 
2010; Pereira-Santaella etal. 2011; Mineo etal. 2012). For 
passive early-type galaxies. X-ray emission is dominated by 
low-mass X-ray binaries (LMXBs) and hot X-ray emitting gas 
(e.g., O'Suflivan etal. 2001; Gilfanov 2004; Boroson etal. 
2011). In this section, we study the normal-galaxy number 
counts as a function of redshift and morphology. 

We measured the number counts of normal galaxies in 
two redshift bins divided at the median redshift Zmedian ~ 0.6. 
In Figure 13, we show the HB and SB number counts 
for normal galaxies at z ^ 0.6 and z ^ 0.6. In the SB, 
we find that low-redshift (z ^ 0.6) normal galaxies domi- 
nate the bright-end counts; however, below flux levels of 
ss(l-2) xlO"'^ ergs cm"^ s"', high-redshift (z ^ 0.6) nor- 
mal galaxies become the dominant galaxy population. By 
contrast, normal-galaxy number counts in the HB are dom- 
inated by low-redshift normal galaxies across the full flux 
range, with ^ 20% of the number counts being from the high- 
redshift population. 

To characterize the number counts of normal galaxies di- 
vided by morphology, we made use of the multiwavelength 
data to classify the normal galaxies broadly as either late-type 
or early-type galaxies. We matched the 169 X-ray detected 
normal galaxies to the optical source catalogs constructed by 
Xue etal. (2010), which provide rest-frame optical magni- 
tudes based on SED fitting. To distinguish between relatively 
blue late-type galaxies and red early-type galaxies, we made 
use of the rest-frame U -V color As described by Bell et al. 
(2004), the U and V bandpass pair straddle the 4000 A break 
and as a color provide a first-order indicator of the stellar age 
of the galactic stellar population. We applied the empirically 
calibrated redshift-dependent color division calculated in § 5 
of Bell et al. (2004) to separate active blue and passive red 
galaxy populations. This division was set at 

(t/-V)rest= l-15-0.31z-0.08(My-l-20.7), (8) 
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Fig. 1 1. — Similar to Figure 9, but for AGNs divided into ranges of intrinsic column density (log[A'H/cm ^] < 22 [open red circles], logiNyi/cm ^] = 22- 
23 [open green triangles], log[A'H/cm"^] > 23 [open blue squares], and no measured Nw due to lack of redshift measurements [filled gray circles]). As in 
Figure 9, AGN number-count predictions from G07 have been provided as dotted and solid curves, corresponding respectively to the inclusion and exclusion of 
an exponential decline in the AGN space density at z > 2.7 (model curve colors correspond to symbol colors). [A machine-readable table of the Nn-divided AGN 
cumulative number counts data for all four bandpasses is provided in the electronic edition.] 




10"" 10"'" 10"'^ 10"'" 10"'^ 10"'° 10"'° 10"'" 10" 
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Fig. 12. — Same as Figure 9, but for AGNs divided into ranges of intrinsic 0.5-8 keV luminosity (log[Lx/ergs s '] < 42 [open red circles], 
log[Lx/ergs s"'] = 42-A3 [open green triangles], log[Lx/ergs s"'] = 43^4 [open blue squares], log[Lx/ergs s"'] > 44 [open orange stars], and no measured 
Lx due to lack of redshift measurements [filled gray circles^). As in Figure 9, AGN number-count predictions from GO? have been provided as dotted and solid 
curves, corresponding respectively to the inclusion and exclusion of an exponential decline in the AGN space density at z > 2.7 (model curve colors correspond to 
symbol colors). [A machine-readable table of the Lx-divided AGN cumulative number counts data for all four bandpasses is provided in the electronic edition.] 
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such that galaxies redward and blueward of this division are 
candidate early-type and late-type galaxies, respectively. We 
note that a single color division will not perfectly isolate the 
early-type and late-type galaxy populations. For example, 
edge-on or dusty late-type galaxies may be artificially red- 
dened and therefore classified as being early-type galaxies 
(see, e.g., Cardamone etal. 2010; Lusso etal. 2011). Fur- 
thermore, intrinsic scatter in the optical color distributions 
near the color cut will sometimes displace galaxies to the 
wrong side of the color division (e.g., "green" early-type 
galaxies that undergo low-levels of recent star-formation). We 
therefore visually inspected HST z^so-band images of all 169 
galaxies to identify obvious cases where rest-frame colors 
provided incorrect morphological classifications. We found 
that 25 galaxies («15%) had obviously incorrect initial mor- 
phological classifications that were corrected. 

Out of the 169 normal galaxies, we classified 123 as late- 
type and 43 as early-type galaxies. In Figure 14, we show 
example HST three-color (B435, Vfiof,, and z^so) images of 10 
early-type {top) and 10 late-type {bottom) galaxies that were 
in the normal galaxy sample. Following the prescriptions 
discussed in § 2, we computed number counts for each of 
these populations. In Figure 15, we show the number counts 
of the normal galaxies broken down by optical morphology. 
The normal-galaxy number counts for early-type and late- 
type galaxies are comparable at the brightest SB and HB flux 
levels. Progressing to lower fluxes, however, we find that the 
late-type galaxy population quickly rises and dominates the 
normal-galaxy number counts, and at the survey flux limits, 
late-type galaxies comprise Ri80% and Ri90% of the normal 
galaxy number counts for the SB and HB, respectively (see 
lower panels of Fig. 15). 

3.4. Contributions to the Cosmic X-ray Background 

With the number-count estimates derived above, we 
can measure the corresponding contributions each source 
type and subtype make to the extragalactic CXRB. Here- 
after, we adopt CXRB intensities JIcxrb of (8.15 ± 
0.58) X 10-12 ergs cm'^ s"' deg-^ and (1.73 ± 0.23) x 
10~'i ergs cm-2 s"^ deg"^ for the SB and HB, respectively. 
These CXRB intensities were computed by summing com- 
ponents from (1) the wl Ms CDF-S unresolved background 
intensity (ilunres; Hickox & Markevitch 2006), (2) the in- 
tensity of faint «1 Ms CDF-S sources (fifaint; below ~5 x 
10"'^ ergs cm~2 s~' and «1.4 x 10"^^ ergs cm"^ s"^ for the 
SB and HB, resepectively) as derived by Hickox & Marke- 
vitch (2006; see their Table 5), and (3) the bright-source in- 
tensity (ilbiight) derived from the best-fit number-counts re- 
lations by Kim etal. (2007), which are based on ~5500 
X-ray sources from the Chandra Multiwavelength Project 
(ChaMP; e.g.. Green et al. 2004; Kim et al. 2004a,b). In the 

bright-source case, f^bright = S'{dN/dS')dS', where 5fai„, 
is the brightest flux used to compute the faint-end counts, 
^bright = 10"'' ergs cm"2 s"' (for both the SB and HB), and 
dN/dS' is the best-fit differential number counts model from 
Table 3 of Kim etal. (2007). Thus, we derived ricxRB = 
flumes + ^faint + ^^bright for the SB and HB . For the FB and UHB , 
we adopted values of (2.54 ±0.24) x 10~" ergs cm"^ s"' deg"^ 
(the sum of the SB and HB CXRB intensities) and (1.04 ± 
0.23) X 10~" ergs cm"^ s~' deg"^, respectively. The UHB 
CXRB intensity was computed by converting the HB inten- 
sity to UHB assuming a power- law SED with F = 1.4 (see 
Morettietal. 2009). 



In Table 2 (col. 1 1-14), we provide the intensities and frac- 
tional CXRB contributions that sources in various subcate- 
gories in the «4 Ms CDF-S provide. We find that the vast 
majority of the CXRB in each band is expected to be produced 
by AGNs, with normal galaxies producing only a small frac- 
tion (r:!0.5-3.0%) of the intensities. These results are similar 
to those found by B04, although revised to fainter flux levels. 

Due to the relatively small solid-angle coverage of the 
CDF-S survey («465 arcmin^), we are unable to characterize 
the number counts for the relatively rare X-ray bright source 
population (i.e., S ^ [1-10] x lO"'"* ergs cm"-^ s~'). To char- 
acterize the bright-end counts for each source type and sub- 
type listed in Table 2 (col. 2) would require a wider survey 
with equivalent multiwavelength data to that used here for 
the «4 Ms CDF-S. Despite this Umitation, we can utilize 
our number-counts estimates and the total bright-end CXRB 
intensities (i.e., r^bnght) to estimate the total resolved CXRB 
intensities. In the first row of Table 2, we provide the bright- 
end corrected resolved source intensities and fractional con- 
tributions to flie CXRB. We find that «76-88% of the CXRB 
intensity can be attributed to X-ray point sources, with the re- 
solved fraction appears to increase with median bandpass en- 
ergy; however, a constant or decreasing resolved fraction with 
energy is not formally ruled out. Comparisons with past stud- 
ies indicate a variety of levels of consistency with this result. 
For example, the studies of Hickox & Markevitch (2006) and 
Kim et al. (2007), which were used to estimate the total CXRB 
intensities and bright-end corrections for the CDF-S, find sim- 
ilar resolved CXRB fractions (using « 1-2 Ms depth CDF 
data) over the same energy ranges studied here. However, the 
studies of Morretti etal. (2003), B04, Worsley etal. (2005), 
and G08 find that the resolved CXRB fraction appears to de- 
cUne somewhat with energy, using data that reaches similar 
Chandra depths («l-2 Ms). These studies adopt results from 
either Morretti et al. (2003) and/or De Luca & Molendi (2004) 
when estimating the CXRB intensity and bright-end counts. 
We note, however, that even in these studies, the apparent 
decline in resolved CXRB with energy is less significant in 
the CDF-S itself (see, e.g.. Fig. 2 of Worsley etal. 2005). 
We estimate that between the wl-2 Ms and ^4 Ms depths, 
the fraction of the resolved CXRB has increased by only 
« 1-2%, which is much smaller than the error bars of our mea- 
surements. Therefore, differences between previous resolved 
CXRB fractions are primarily related to either differing as- 
sumptions about the CXRB intensity and bright-end counts 
or field-to-field variations. In support of the latter point, LOS 
find that the CDF-S number counts are lower than those of the 
CDF-N by k.25% at the faintest fluxes. 

4. COMPARISONS WITH MODELS 

4.1. The AGN Number Counts and the Evolution of Accreting 

SMBHs 

The redshift-dependent AGN number counts (Figs. 9 and 
10) provide a directly observable signature of the evolution of 
AGN activity in the Universe. Since AGNs are the brightest 
X-ray detected population, many previous investigations have 
focused on measuring their cosmic evolution (see Brandt & 
Hasinger 2005 for a review), and in particular, the evolution of 
their X-ray luminosity function (XLF; e.g., Cowie et al. 2003; 
Ueda et al. 2003; Hasinger et al. 2005; Silverman et al. 2008; 
Ebrero etal. 2009; Aird etal. 2010). For example, Hasinger 
etal. (2005) made use of a variety of ROSAT, XMM-Newton, 
and Chandra surveys with highly-complete optical/near-IR 
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redshift measurements to compute the observed evolution 
of the AGN XLF. This investigation found that X-ray se- 
lected AGNs appear to undergo luminosity-dependent den- 
sity evolution in their populations, where the space density 
of low-luminosity (Seyfert-type) AGNs peaks at z ^ 1 and 
the highly luminous (QSO-type) AGN space density peaks 
at z ~ 2. These direct measurements of the observed XLF 
evolution of AGN activity in the Universe provide a useful 
constraint on the SMBH accretion history; however, the ob- 
served X-ray bandpass used to make these measurements (i.e., 
«0.5-2 keV) is susceptible to extinction and therefore popu- 
lations of AGNs with intrinsic X-ray absorption columns of 
A^H ^ 10^^ cm"^ will have under-represented intrinsic accre- 
tion activity. 

Using the unresolved hard X-ray (?s3-100 keV) back- 
ground intensity spectrum as a constraint on the obscured 
(A^H ~ 1022_io24 cm-2) and Compton-thick (A^h ^ lO^'* cm'^) 
AGN populations, Gilli etal. (2007; hereafter, GOV) used the 
observed X-ray luminosity function evolution from Hasinger 
et al. (2005), local distributions of the luminosity-dependent 
obscured-to-unobscured AGN fraction, and a model distri- 
bution of AGN population column densities to form phe- 
nomenological models describing the intrinsic evolution of 
the AGN X-ray luminosity function (see also Ueda et al. 2003; 
Draper & Ballantyne 2009; Treister etal. 2009 for additional 
prescriptions). These models make direct predictions for the 
number count distributions of sources with different ranges 
of z, Nu, and intrinsic 0.5-8 keV luminosity LxC3 and can 
be directly compared with our observed AGN number-count 
subpopulations (see Figs. 9-12). 

In Figures 9-12, we have highlighted the number-count 
predictions for z, A^h, and Lx selected AGNs by G07 appro- 
priate for the Hasinger et al. (2005) XLF as observed (solid 
curves) and after applying a high-redshift (z > 2.7) exponen- 
tial decline to the XLF [dotted curves; see Schmidt et al. 1995 
for motivation). In general, we found better agreement by not 
including the exponential declining term (i.e., the solid curves 
in Figs. 9-12). Although, since the majority of the AGNs in 
this study are at z ^ 2.7, there is generally little difference be- 
tween the models that include and do not include an exponen- 
tially declining term. For clarity in comparisons, in Figs. 9, 
11, and 12 we have provided lower panels showing the ra- 
tio between our observed number counts and the G07 models 
without the exponential declining term. 

For redshift-selected number counts, we found reasonable 
agreement (within a factor of w2) between the observed num- 
ber counts and the G07 model across the entire SB and HB 
flux ranges (see Fig. 9). We note that a non-negligible frac- 
tion of the AGN population does not have redshifts available. 
Since the multiwavelength coverage over the entire CDF-S is 
extensive and deep, sources without redshifts (neither spec- 
troscopic nor photometric) have very faint optical/near-lR 
counterparts, and are are therefore good candidates for high- 
redshift (z > 3) AGNs. If we add the number counts from 
AGNs without redshifts to the z > 3 number counts, we 
find improved agreement with the G07 model prediction (see 
Fig. 10 diamonds). For comparison, we have also plotted the 
predictions from the recent work by Aird et al. (2010), which 
is based on 2-10 keV XLF measurements from the CDFs and 
the AEGIS -X surveys. The Aird etal. (2010) predictions are 
also consistent with our z > 3 AGN number counts; however. 

The G07 phenomenological model number-counts predictions can be 
found at http://www.bo.astro.it/~gilIi/counts.html 



if all of the AGNs with unknown redshifts are in the z > 3 
range, our number counts become more consistent with the 
G07 model. Regardless of the redshifts of the sources lacking 
redshift identifications, our number counts appear to be a fac- 
tor of ^ 2 higher than those predicted by the G07 model with 
exponential decline. 

For the number counts predictions based on A^h, we find 
significant discrepancies (factors of sal. 5-3) between those 
observed and those predicted by the G07 model (see bot- 
tom panels of Fig. 11). This situation is similar regardless of 
whether the exponential decline is applied or not. We suspect 
that some of the disagreement here will be due to (1) cosmic 
variance of the populations, (2) our simplistic calculation of 
A^H based on band ratios (see discussion in §3.2), which differs 
from the spectral models from G07 (see their Fig. 1), (3) the 
distribution of A^h values for sources without redshifts, and 
(4) imperfections in the G07 model descriptions. Regarding 
point (3), we note that sources with unknown redshifts have 
a median effective photon index of Fmedian ~ 0.7, which is 
somewhat harder (softer) than the Fmedian ~ 1 .0 (Fmedian ~ 0.4) 
for sources with log(A^H/cm"2) = 22-23 (log[A^H/cm"2] > 23). 
Therefore, we expect that most of the sources with unknown 
redshifts will have log(A^H/cm"^) ^ 23 regardless of their red- 
shifts. This would improve the agreement between observed 
and model number counts for log(A^H/cm"^) ^ 23. Despite 
the apparent disagreements and limitations, we do find basic 
agreement in the general trend: more heavily obscured AGNs 
have larger contributions at fainter fluxes. 

For number counts divided by intrinsic 0.5-8 keV luminos- 
ity Lx, we find good (within a factor of ?s2) basic agreement 
between our measurements and the G07 model (see Fig. 12); 
however, this is not true for the decline model (see dotted 
curves in upper panels of Fig. 12). 

In summary, we find that the G07 phenomenological mod- 
els provide a good description of the redshift, column density, 
and intrinsic luminosity distributions to the extent that we can 
constrain these values observationally. Better measurements 
of the column densities and Lx would be needed to improve 
these models further However, unless a significant fraction 
of objects with photometric redshifts z > 3 prove to be lower 
redshift contaminants, we do find that applying an exponential 
decline to the AGN XLFs at z > 2.7, under the G07 param- 
eterizations, provides a poor characterization of the redshift 
and Lx divided number counts. 

4.2. The Rapidly Increasing Normal Galaxy Counts 

In § 3.3, we presented number-counts estimates for the 
normal-galaxy population and subpopulations separated by 
redshift and galaxy morphology. Given that these X-ray 
source populations are uniquely accessible to the ultradeep 
survey regime probed by the Chandra Deep Fields, the evo- 
lution of the normal-galaxy XLFs have not been as well char- 
acterized as those of AGN populations (see § 4.1). However, 
initial investigations have placed first-order constraints on the 
evolution of normal-galaxy XLFs out to z ~ 0.75 (e.g., Nor- 
man etal. 2004; Georgakakis etal. 2007; Ptak etal. 2007; 
Tzanavaris & Georgantopolous 2008). These investigations 
have provided evidence for rapidly increasing late-type galaxy 
0.5-2 keV XLF that is consistent with pure luminosity evolu- 
tion, where L* cx (l+z)*", and w 1.5-3. By contrast, the 
early-type galaxy 0.5-2 keV XLF appears to be changing 
more slowly with redshift, and in the Tzanavaris & Georgan- 
topolous (2008) study, the XLF is consistent with being con- 
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Fig. 13. — Top Panels: Cumulative nomial-galaxy number counts (filled black circles) in the SB (left) and HB (right) broken down by contributions from 
galaxies at redshifts z ^ 0.6 (open orange stars) and z 0.6 (open green upside-down triangles). The redshift division between the two normal-galaxy samples 
was made at the median galaxy redshift. The dotted curves in the left panel show the predicted z ~ 0.6 (orange dotted curx'e) and z 0.6 (green dotted curve) 
number counts based on the observed XLF pure-luminosity evolution parameterization derived by Ptak etal. (2007) using observed XLFs out to z 0.75. 
At z 0.6, these predictions are based on an extrapolation of the XLF parameterization out to z~2. The solid curves in the left and right panels show the 
predicted number counts for z ^ 6 (orange solid cur\'es) and z ^ 0.6 (green solid cun'es) galaxy populations based on the observed evolution of stellar mass 
functions converted to XLFs using X-ray scaling relations (see §4.2 for details). Middle Panels: The fractional contributions that z ~ 0.6 and z ^ 0.6 galaxy 
populations make to the total galaxy number counts. Bottom Panels: The ratio between the data and predictions based on the solid curves in the top panels. [A 
machine-readable table of the redshift-divided normal galaxy cumulative number counts data for all four bandpasses is provided in the electronic edition.} 



stant or decreasing with redshift (the value of b ranges from 
« -2.3 to +2.6 in these studies). However, these results are 
based on shallower data than the ?s4 Ms CDF-S, and thus far 
provide only estimates for how the bright-end of the galaxy 
XLFs evolved out to z ss 0.75. In the «4 Ms CDF-S, w29% 
of the normal galaxies are at z ^ 0.75, and the majority of 
these sources are only accessible by this survey. In the para- 
graphs below, we interpret the number counts of the normal 
galaxy subpopulations (e.g., see Figs. 13 and 15) in the con- 
text of what we expect from previous estimates of the galaxy 
XLF evolution and the evolution of galaxy properties (e.g., 
SFR and stellar mass) and X-ray scaling relations (e.g., the 
X-ray/SFR relation). 

Following the formalism of Ranalli etal. (2005), the 
normal-galaxy number counts can be expressed as the follow- 
ing double integral: 

(9) 

where K^^^ = 3.05 x 10""* converts from sr~' to deg~^, ipx is 
the redshift-dependent XLF, -^^^ is the cosmology-dependent 
differential volume element (comoving volume per unit red- 
shift per unit solid angle), and [Zmm, Zmax] and [L,nin(5), 
Lmax] are respectively redshift and luminosity integration lim- 
its. The value of Lmin(5') is the rest-frame luminosity corre- 
sponding to a flux S and redshift z- In the upper-left pan- 
els of Figures 13 and 15, the dotted curves show 0.5-2 keV 
number count predictions based on previous measurements 



of the evolution of galaxy XLFs for redshift and morphol- 
ogy divided galaxy populations, respectively. These curves 
were computed using the Ptak et al. (2007) redshift-dependent 
Schechter parameterizations of the observed galaxy XLFs, 
which assume pure-luminosity evolution of the form L* cx 
(1 +z)^'^ and (1 +z)''* for late-types and early-types, respec- 
tively. We also experimented with using a log-normal param- 
eterization of the galaxy XLF and found similar results. When 
estimating galaxy number counts predictions, we used equa- 
tion 9 and extrapolated the Ptak et al. (2007) parameteriza- 
tions out to z « 2. The Ptak et al. (2007) XLFs were directly 
measured for galaxies out to z « 0.75 (using galaxies with red- 
shifts as large as z ~ 1 .2) and are in good agreement with other 
galaxy XLF measurements from the literature (e.g., Norman 
etal. 2004; Tzanavaris & Georgantopolous 2008). The un- 
certainties on the normalization of the XLF parameterization 
is a factor of ~2. We find that, within the uncertainties, the 
XLFs provide reasonable predictions for the z ^ 0.6 (Fig. 13, 
orange dotted curve and stars) and late-type galaxy (Fig. 15, 
blue dotted curve and squares) number counts. For these two 
samples, the largest differences come from the z ^ 0.6 sam- 
ple at 0.5-2 keV fluxes ^3 x 10"'^ ergs cm"^ s"', which 
may be due to differences in galaxy classification of X-ray 
detected sources, cosmic variance, and statistical errors as- 
sociated with the XLF fits. We find that the extrapolated 
XLFs appear to overpredict the z^ 0.6 (Fig. 13, green dot- 
ted curve and upside-down triangles) and early-type galaxy 
(Fig. 15, red dotted curves and triangles) number counts by 
more than a factor of two, suggesting that our extrapolation of 
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Fig. 14. — Thi'ee-color (B435, V^^, and zgso) images of a sample of X-ray detected nonnal galaxies ordered by redshift. The top panels show 10 example 
late-type galaxies and the bottom panels show 10 example early-type galaxies. The angular extent of each image has been chosen to subtend a physical distance 
of Ri30 kpc on a side in the rest-frame of each galaxy. The redshift of each source is indicated in the upper left and a 1 arcsec bar has been placed in the lower-left 
hand comer to show angular scale. The images are centered on the coordinates corresponding to each galaxy, which were taken from columns 18 and 19 of 
Table 3 in XI 1. These positions were obtained primarily from VLT WFI observations, and therefore small ( ^ 0.1-0.2 arcsec) offsets are sometimes seen in these 
HST images (note the angular scale). 



the early-type galaxy XLF evolution out to z w 2 is an over- 
estimate. Therefore, it is possible that the rate by which the 
early-type galaxy XLF rises with redshift decreases somewhat 
going from z ~ 0.75-2 and/or the adopted rate of increase is 
too high (e.g., see Tzanavaris & Georgantopolous 2008 who 
find a lower rate of increase). To obtain independent predic- 
tions of the galaxy XLFs and number counts for galaxies out 
to z w 2, we can use the combination of X-ray scaling rela- 
tions (e.g., the X-ray/SFR relation) and the redshift evolution 
of galaxy physical properties (e.g., SFR and stellar mass) that 
are measured using multiwavelength observations. 

Recent X-ray stacking investigations of large normal- 
galaxy populations selected by morphology and galaxy phys- 
ical properties (e.g., SFR and stellar mass) have now provided 
insight into how the X-ray emission from normal galaxies 
evolves relative to galaxy physical properties out to z « 1- 
4. Lehmer et al. (2008) reported that late-type star-forming 
galaxies, the most numerous galaxy population in the Uni- 



verse, have mean Lx/SFR values that are roughly constant out 
to z « 1 .4. Further X-ray stacking of distant z ~ 1 .5^ Lyman 
break galaxies appears to show similar mean Lx/SFR values 
(e.g., Brandt etal. 2001; Lehmer etal. 2005; Laird etal. 2006; 
Cowie etal. 2011; Basu-Zych etal. 2012, in-prep). Similarly, 
Vattakunnel etal. (2012) have studied z ^ 1.2 normal galax- 
ies in the «4 Ms CDF-S at X-ray and radio wavelengths, and 
found that the Lx-i'i.4GHz correlation for these galaxies is 
consistent with that expected from the local Lx/SFR relation, 
albeit with significant scatter (wO.4 dex). Additional indirect 
support for a roughly constant Lx/SFR ratio with redshift has 
been provided by Dijkstra et al. (201 1) who used the observed 
SFR density evolution of the Universe to show that the unre- 
solved SB CXRB can be fully explained by a Lx/SFR ratio 
oc (1 +z)*, where b is constrained to be less than 1.4. 

For early-type galaxies, Lehmer et al. (2007) and Danielson 
etal. (2012) showed that the soft-band X-ray luminosity per 
unit Z?-band luminosity (Lx/Lb) appears to undergo only mild 
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properties to make predictions for the evolution of the normal- 
galaxy XLFs and number counts for galaxies out to z w 2. 

Following Avni & Tananbaum (1986), Lpx can be estimated 
through the following transformation: 

/•oo 

(px(logLx,z)= / (^i)(logi?,z)f(logLx|logi9)t/logz?, 



Fig. 15. — Top Panels: Similar to Figure 13, but for normal galaxies divided into late-type {open blue squares) and early-type {open red triangles) galaxy 
populations. The classification of galaxy morphology is explained in §3.3. The dotted curves in the upper left panel show the predicted late-type (blue dotted 
ciin'e) and early-type {red dotted curve) galaxy number counts based on the observed XLF pure luminosity evolution parameterization derived by Ptak etal. 
(2007) using observed XLFs out to z « 0.75. These predictions are based on extrapolation of the XLF parameterization out to 7,~ 2. The solid curves in 
the left and right panels show the predicted number counts for late-type {blue solid curiae) and early-type {red solid curve) galaxy populations based on the 
observed mass functions converted to XLFs using X-ray scaling relations (see §4.2 for details). Middle Panels: The fractional contributions that late-type and 
early-type galaxies make to the total galaxy population. Bottom Panels: The ratio between the data and predictions based on the solid curves in the top panels. 
The reasonable agreement between the number-counts data and models, both here and in Figure 13, indicates that the rapidly rising galaxy counts are largely 
attributed to the evolution of galaxy physical properties (e.g., stellar mass and star-formation rate) with little redshift evolution in scaling relations between these 
physical properties and X-ray power output. [A machine-readable table of the morphology-divided normal galaxy cumulative number counts data for all four 
bandpasses is provided in the electronic edition.} 

evolution (cx [1 -l-z] to z w 1 .2 for the most optically lumi- 
nous galaxies, which is expected to be due to X-ray emitting 
hot gas. In contrast, however, lower-luminosity early-type 
galaxies have been found to have rising mean Lx/Lb values 
with increasing redshift, which is expected to be due to a fad- 
ing LMXB population toward the present day (Lehmer etal. 
2007; Homschemeier etal. 2012, in prep). Similar expecta- 
tions have been found from X-ray binary population-synthesis 
models, which predict an evolution of Lx/Lb (x (1-l-z)^ at 
1 (Fragos etal. 2012, in prep; see also Ghosh & White 
2001). 

The above investigations provide estimates for how scal- 
ing relations between X-ray luminosity and galaxy physical 
properties have evolved out to z « 1-4, where the majority 
of the w4 Ms CDF-S normal galaxies are detected. It seems 
that only the scaling relation between LMXB emission and 
stellar mass undergoes significant evolution with cosmic time; 
however, these populations are expected to dominate only in 
early-type galaxies, which make up a minority of the X-ray 
detected sources in the w4 Ms CDF-S. It is therefore likely 
that the evolving normal-galaxy XLF and rising galaxy num- 
ber counts, observed here and in previous studies, is largely 
due to the rapid evolution of galaxy physical properties. Re- 
cent multiwavelength surveys (e.g., COSMOS; Scoville etal. 
2007) have provided significant new constraints on the evo- 
lution of optical/near-lR luminosity, SFR, and stellar-mass 
functions. These measurements can be combined with scal- 
ing relations between X-ray luminosity and galaxy physical 



(10) 

where -d represents a galaxy physical property (e.g., optical 
luminosity, stellar mass, or SFR) by which ip^ilog-d^z) has 
been measured, and P{logLx\log-d) is the probability distri- 
bution for observing Lx given the value -d. 

We make use of the redshift-dependent stellar-mass func- 
tions measured by llbert etal. (2010) for late-type and early- 
type galaxies out to z ~ 2 (i.e., d = Mi,). To convert 
(^M,(logMj,,z) to </Jx(logix,z) via equation 9 requires knowl- 
edge of the transformation between X-ray luminosity and stel- 
lar mass (i.e., P[logLx| logM,t]). This function is expected to 
depend on galaxy morphology. In its most general form, we 
can express the transformation function as follows: 



P(logLx|logMJ = 



1 



27r(T 



exp 



(logLx,mod(M*)-logLx) 
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(11) 

where Lx.modiM*) is the predicted X-ray luminosity given a 
value of the stellar mass and cr represents the scatter in the 
relation. For normal late-type galaxies in the local Universe, 
the X-ray luminosity has been shown to correlate strongly 
with SFR (e.g., Persic & Rephaeli 2007; Lehmer etal. 2010; 
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Pereira-Santaella etal. 2011; Symeonidis etal. 2011; Mineo 
etal. 2012), and as discussed above, this correlation is ob- 
served to hold out to z ^ 1. We made use of the Lehmer etal. 
(2010) X-ray/SFR relation for late-type star-forming galaxies 
in the local Universe: 

^2-10 keV = aM* + /3SFR (12) 

where a = 9.08 x lO^^ ergs s"' and ;3 = 1.62 x 

10^^ ergs s"' (M0yr"')"', Mi, is in units of Mq, SFR is in 
units of Mq yr"', and -L2-10 keV is in units of ergs s"'. Stud- 
ies of local and distant galaxy populations have shown that 
the SFR for late- type star-forming galaxies is further strongly 
correlated with M^, out to z « 2 (e.g., Daddi etal. 2007; Elbaz 
etal. 2007; Peng etal. 2010), albeit with a rapidly evolving 
normalization (i.e., mean SFR/M*) with redshift (e.g., Zheng 
et al. 2007; Karim et al. 201 1; Cen 201 1). To convert SFR to 
for late-type star-forming galaxies, we adopted the rela- 
tion provided by Karim et al. (201 1): 

SFR/M, = co(l +z)^-^(M,,/10"Mo)*™ (13) 

where cq « 0.0263 Gyr"' and /3sfr ~ -0.4 are fitting con- 
stants. Combining this relation with equation 12, we arrive 
at the following expression: 

L2-10 kev = aM* + 7(l +zf-^M^J (14) 

where 7 = 1.07 x 10^-^, M^, is in units of Mq, and L2-10 keV is 
in units of ergs s~^ When computing number-count predic- 
tions, we converted the 2-10 keV luminosities provided here 
to other bandpasses assuming a power-law SED with F = 1 .9, 
the average stacked SED of the X-ray detected normal galax- 
ies (Young etal. 2012). The overall scatter (accounting for 
both SFR/M,t and Lx/SFR relations) in this relation is esti- 
mated to be (T « 0.4 dex. 

For normal early-type galaxies, the X-ray luminosity has 
been shown to correlate with /T-band luminosity, Lk, which 
provides a direct proxy for galaxy stellar mass (e.g., Gil- 
fanov 2004; Boroson etal. 2011). We converted the Boro- 
son etal. (2011) Lx/Lk relations for early-type galaxies 
to LxIMi, relations assuming a single mass-to-hght ratio 
log /La: « -0.2 (where Mi, and Lk are in solar units) char- 
acteristic of early-type galaxies. These relations are separated 
into contributions from LMXBs and hot X-ray emitting gas: 

logLo.3-8 kev(LMXB) w 29.2+logM„-h21og(l +z) 

logZ-o.3-8 kev(gas) « -6.07+4.03 logM^, (15) 

where the 21og(l+z) term accounts for the expected red- 
shift evolution for LMXBs in early-type galaxies (see dis- 
cussion above). When computing number-counts predictions, 
we converted the 0.3-8 keV relations to other bandpasses as- 
suming a power-law SED with F = 1.8 for LMXBs and an 
apec plasma model with an Lo.3-8 keV dependent temperature 
{Jx cx ^03!^ i^gv) solar abundances for the hot gas. The 
scatter in the LMXB and hot-gas relations are estimated to 
be «0.2 dex and «1.0 dex, respectively (see Boroson etal. 
2011). 

The relations expressed in equations 14 and 15, combined 
with the stellar-mass functions from Ilbert et al. (2010), were 
folded into equation 9 to obtain redshift-dependent XLFs for 
late-type and early-type galaxies, respectively. We measured 
the number counts of both populations in two redshift bins 
by integrating equation 9 over the redshift intervals z ^ 0.6 
and z ^ 0.6 and summing the contributions from each galaxy 



morphology type. In Figures 13 and 15, we show our model 
number counts in bins of redshift and morphology, respec- 
tively, and in the bottom panels, we plot the ratio of our data 
to the model. We remind the reader that the models im- 
plemented here are not fits to the data in any way. For the 
redshift-divided normal galaxy samples, we find good agree- 
ment between our data and the models at the SB and HB flux 
limits; however, for the number counts for normal galaxies 
at z ^ 0.6, we find significant differences at the bright-end 
(/o.5-2keV and /i-gkev ~ 10"'^ ergs cm'^ s"'). The bright- 
end number counts for z ^ 0.6 galaxies are based on only 
wl-5 objects; therefore, the differences with the models may 
be due to either cosmic variance or misclassification of some 
fraction of these sources (i.e., if some of these normal galax- 
ies are actually AGNs). Furthermore, the general trends seen 
in the data and models (e.g., the increasing contributions and 
faint-end dominance in the SB of z ^ 0.6 galaxies) appear to 
be broadly consistent with each other. For the morphology- 
selected number counts, we find very good agreement (less 
than a factor of «1.5 difference) between the data and the 
model with the exception of the HB early-type galaxy counts. 
We note, however, that the HB number counts measurements 
for early-type galaxies are based on small numbers of galax- 
ies (four galaxies) and are sensitive to the classifications dis- 
cussed in §3.1 above. 

In summary, we find that our models provide good char- 
acterizations of the X-ray number counts in redshift bins for 
both late-type and early-type galaxy populations down to the 
flux limits of our survey (with the exception of z 0.6 nor- 
mal galaxies with bright SB and HB fluxes and HB early-type 
galaxy counts; however, see discussion above). For late-type 
galaxies, this result suggests that X-ray scahng relations in 
the local universe appear to hold out to at least z « 1-2, in 
agreement with previous investigations. For early-type galax- 
ies, these results are consistent with there being little redshift 
evolution in how hot gas X-ray emission scales with stellar 
mass and modest redshift evolution (oc [1 + z]^) in the LMXB 
luminosity per stellar mass. 

5. THE IMPENDING DOMINANCE OF NORMAL GALAXIES: 
PREDICTIONS FOR FUTURE DEEPER X-RAY OBSERVATIONS 

The number-counts results presented above show that nor- 
mal galaxies play an increasingly dominant role in the extra- 
galactic X-ray source population going to fainter fluxes. Our 
observations show that the cumulative SB and FB number 
counts of normal galaxies are on the brink of surpassing those 
of AGNs (see Fig. 5), and the slope of the galaxy number 
counts (dN/dS) already exceeds that of AGNs (see Fig. 6). 
At the SB flux limit of the ^4 Ms CDF-S normal galaxies are 
reaching the unprecedented source density of « 12,700 deg~^, 
which constitutes «46% of the total number counts. Given the 
comparable source densities of normal galaxies and AGNs at 
the survey flux limit and the relatively sharp faint-end rise in 
the galaxy number counts versus AGNs (A'[> 5] power-law 
slope of /3s^'- 1 « 1.2 versus /3f™- 1 w 0.5; see Table 1), it 
is clear that going to fainter Chandra fluxes beyond w4 Ms 
will yield larger gains in X-ray source densities than previous 
increases in depth. XI 1 showed that exposures of the CDF-S 
region will yield substantial gains in sensitivity. For example, 
afeasible «10Ms exposure would permit a factor of «1.8 in- 
crease in depth (in terms of flux limit) over the current ~4 Ms 
exposure. Therefore, exploring the normal-galaxy dominated 
regime of the number-counts distribution is a sensible goal 
that may be pursued in future surveys. At present, the only 
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Fig. 16. — Measured SB number counts (filled black circles) with con- 
tributions from AGNs {open blue triangles) and normal galaxies (open red 
squares) shown. Our best-fit parameterizations from equation 5 for all 
sources, AGNs, and normal galaxies have been shown as solid black, dashed 
blue, and dotted red curves, respectively. An extrapolation of these models 
indicates that the normal-galaxy number counts will overtake those of AGNs 
at SB flux levels of S:i(3-5) X 10""* ergs cm"^ s"' . The ultimate flux lim- 
its of the Si4 Ms CDF-S, a potential Ms CDF-S, and future fa4 Ms 
SMART-X exposures have been indicated. These extrapolations indicate that 
a fa 10 Ms CDF-S would have SB number counts at the flux limit that reach 
source densities of ai46,000 deg"^ and are dominated by normal galaxies. 
More ambitiously, a ft!4 Ms SMART-X exposure would reach source den- 
sities of f»500,000 deg^ and would have an Si90% contribution from nor- 
mal galaxies. For reference, we have plotted the source densities of galaxies 
with i?-band magnitudes brighter than R^2A and R^21 (gray short-dashed 
lines). 

X-ray mission capable of reaching these fluxes is Chandra 
(see below), and this will likely be the case for at least the 
next two decades. 

In Figure 16, we provide an expanded view of the SB num- 
ber counts, and their contributions by extragalactic source 
type (AGNs and normal galaxies), in the faint flux regime. 
Our best-fit number-counts parameterizations (based on equa- 
tion 5) have been indicated. Direct extrapolation of these pa- 
rameterizations indicates that the number counts of normal 
galaxies will overtake those of AGNs just below our current 
flux Hmit at «(3-5) x 10""^ ergs cm"^ s"' . With a «10 Ms ex- 
posure in the CDF-S, we could study number counts down to 
a sensitivity limit of «2.8 x 10"'^ ergs cm"^ s"' (see vertical 
dotted line in Fig. 16), although we expect the faintest sources 
will have most-probable fluxes of «5.2 x 10""^ ergs cm"- s"' 
(see §2.1 for distinction). At such a depth, we would ex- 
pect a factor of «1.7 increase in total X-ray source densi- 
ties to «46,000 deg"^ (based on extrapolation of our best-fit 
number-counts parameterization), with normal galaxies mak- 
ing up s=;57% of the limiting source density. After consid- 
ering variations in sensitivity across the field, we expect that 
a 10 Ms exposure would yield wl,020-l,080 total X-ray de- 
tected sources. 

While Chandra is the only operating observatory capable 
of observing to the ultra-deep regime studied here, the re- 



alization of future X-ray mission concepts like Generation- 
X and SMART-X\^ which have much larger light-collecting 
areas and comparable or improved imaging resolution over 
Chandra, would allow for the exploration of new populations 
of extragalactic X-ray sources, including new populations of 
distant (z ^ 0.6) normal galaxies (see, e.g.. Fig. 13). Under its 
current specification, a «4 Ms SMART-X survey would detect 
SB sources with fluxes as low as w3 x 10"'^ ergs cm"^ s"' . Di- 
rect extrapolation of our number counts to these levels predict 
source densities of w500,000 deg"-, with «90% contribution 
from normal galaxy populations. 

6. SUMMARY 

In this paper, we have provided new quantitative measure- 
ments of the X-ray number counts in the w4 Ms CDF-S using 
fourbandpasses: 0.5-2 keV, 2-8 keV, 4-8 keV, and 0.5-8 keV 
(SB, HB, UHB, and FB, respectively). Our analyses focus 
on the faintest flux regimes only accessible by such ultradeep 
Chandra surveys. We draw the following key conclusions: 

• We make use of a Bayesian approach, flux probabil- 
ity distributions, and maximum-likelihood techniques 
to measure number counts for sources detected in the 
w4 Ms CDF-S down to SB, HB, UHB, and FB flux 
limits of 5.1 X 10"l^ 3.7 x 10"'^ 4.6 x 10"'^ and 
2.4 X 10"'^ ergs cm"^ s"', respectively; these are factors 
of sal. 9^. 3 times fainter than previous number counts 
investigations. 

• At the flux limits of our survey, we reach source den- 
sities of 27,800, 10,500, 8,400, and 22,600 deg"^ in 
the SB, HB, UHB, and FB, respectively. AGNs are 
the majority contributors to the total number counts 
at all fluxes in all four bands (especially in the HB 
and UHB). In the SB, we reach AGN sky densities of 
« 14,900 deg"^, the highest reliable AGN sky density 
measured at any wavelength. 

• The normal-galaxy number counts rapidly rise, com- 
pared with AGNs, and at the SB and FB flux limits, nor- 
mal galaxies make up sa46% and w43% of the number 
counts, respectively. At the limiting SB flux, normal 
galaxies reach a sky density of « 12,700 deg"^, which 
is a factor of Ra4.5 times higher than measured in the 
Ki2 Ms surveys. 

• Stated in terms of the intensity on the sky. X-ray de- 
tected sources in the «4 Ms CDF-S can account for 
«76-88% of the CXRB. We estimate that the increased 
exposure from «l-2 Ms to sa4 Ms has resulted in an in- 
crease of sal -2% resolved CXRB. 

• In die SB and HB, AGNs at high redshifts 
(z ^ 1 .5), with moderate-to-significant obscura- 
tion (A'h ^ 10^^ cm"^), and relatively low intrinsic 
0.5-8 keV luminosity (Lx ^ lO'*'' ergs s"') make in- 
creasingly important contributions to the AGN number 
counts going to fainter fluxes, and at the flux limits, 
these sources make up ^ 40% of the overall AGN 
number counts. These trends are broadly consistent 
with those expected from the phenomenological 
models of G07. 

" See lhttp://www.cfa.harvard.edu/hea/genx/ 1 and http://hea- 
www.cfa.harvard.edu/SMARTX/ for more information about Generation-X 
and SMART-X, respectively. 
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• In the SB and HB, late-type star-forming galaxies with 
z ^ 0.6 provide majority contributions to the normal- 
galaxy number counts, with the exception of the SB 
counts below sslO"'^ ergs cm~^ s~\ where 0.6 
star-forming galaxies dominate. These galaxies are ex- 
pected to produce X-ray emission from high and low 
mass X-ray binaries, hot gas, and young stars. By con- 
trast, passive early-type galaxies, with X-ray emission 
produced by low mass X-ray binaries and hot gas, make 
up only small fractions («5-20%) of the SB and HB 
number counts near the survey flux limits. These trends 
are well described by models that assume (1) the X-ray 
power output for late-type galaxies scales with SFR out 
to z « 2 with no evolution in the scaling relation; and 
(2) the X-ray power output in early-type galaxies is due 
to a non-evolving hot gas emitting component and an 
evolving (by [1 +z]^) low mass X-ray binary component 
that both scale with stellar mass. Therefore, the rapidly 
rising normal-galaxy number counts can be attributed 
primarily to the evolution of the physical properties of 
galaxies (i.e., star formation rate and stellar mass) and 
not X-ray scaling relations. 

• Extrapolation of our number counts to lower fluxes 
suggests that an X-ray observation reaching flux lev- 
els just below those probed by the sa4 Ms CDF-S (at 
SB fluxes of «[3-5] xlO~'** ergs cm"^ s"') would re- 
sult in an X-ray sky with normal galaxies dominating 



the source density at the sensitivity limit. Since the 
normal-galaxy number counts are now comparable with 
AGNs at the ss4 Ms flux limit, and are expected to con- 
tinue to rise rapidly at fainter fluxes, more sensitive 
surveys will have relatively large yields in increased 
source densities. We show that a 10 Ms CDF-S would 
yield « 1,020-1, 080 total sources and would result in a 
limiting SB source density of «46,000 deg"^ with nor- 
mal galaxies dominating the number counts (providing 
«57% of the total number counts) at the flux limit. 
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RECOVERY FRACTION CORRECTIONS 

As noted in §2.2, the XI 1 CDF-S source catalog was generated following a two-step approach, which entailed (1) running 
wavdetect to form an initial list of candidate sources and (2) assessing the probability of detection using AE. Our methods 
described in § 2 for computing number counts generally account for biases due to completeness and the Eddington bias provided 
our catalogs are complete to a specific AE source selection probability. Therefore, we need to account for incompletenesses in 
the number counts that result as a consequence of our two-step cataloging approach. In this section, we summarize our approach 
for computing this correction, which has been implemented in equation 5 as the quantity C'. 

We began by generating 200 mock CDF-S images in the four bandpasses (i.e., 800 images), which each contained 700 X-ray 
sources artificially implanted. These images were constructed as follows. Each source was assigned a random right ascension 
and declination within the area of the ?s4 Ms CDF-S as defined by the exposure maps (see § 3 of XI 1) and was given a random 
number of counts between w 1-100 counts. The counts from each source were added to a blank image canvas using the mar 
(version 4.5) ray-tracing code. In this procedure, the PSF of each source was modeled assuming a roll angle and aim point 
that were selected from one of the 54 observations that make up the cumulative w4 Ms CDF-S exposure (see Table 1 of XI 1). 
Selection of the roll angle and aim point for each source was done probabilistically with the probability of selection being directly 
proportional to the exposure time of each of the 54 CDF-S observations. Mock images for each of the four bandpasses were then 
created by adding background counts from the respective cumulative background images (see §7. 1 of XI 1) to the image canvases. 

For each of the 700 sources in each of the 800 images we performed circular aperture photometry using a circular aperture 
with a radius encompassing w90% of the encircled-energy fraction. The number of counts measured for each source was then 
compared with the number of counts needed for a source detection (see discussion of sensitivity maps in §2. 1) to see if it would 
satisfy the criterion adopted in equation 1 for source detection. We then searched each of the 800 images using wavdetect 
at a false-positive probability threshold of 10"^ and made wavdetect source catalogs for each image. For each bandpass, we 
combined all 200 catalogs of 700 sources and computed the count rate cj) and off-axis angle 9 dependent fraction of sources that 
satisfied our binomial probability selection criterion (i.e., equation 1) that were also detected by wavdetect (hereafter, recovery 
fraction /recov = 1/C')- 

In Figure Al, we show the /recov versus (j> for off-axis angle intervals going to w 9 arcmin, where the image is roughly 
radially contiguous. We found that the data for /recov could be represented successfully using the analytic form /recov(0,^') = 
1/(1 + exp[-S(6){(j> - £^(9)}]), where 6{9) and 1^(6') are fitting constants that vary with off-axis angle and bandpass. In § 2.2, we 
made use of this relation in equation 6, where we adopt C' = 1 //recov We note that this correction has the advantage in that it is 
independent of the Bayesian priors discussed in § 2 and does not depend on the shape of the number counts distribution near and 
below the flux limits. 



See |http://space.mit.edu/cxc/niarx/| for marx simulator details 
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Fig. A1. — Fraction of simulated sources satisfying the binomial probability selection criterion of equation 1 that were also detected using wavdetect at 
a false-positive probablity threshold of 10"^ (i.e., f^ccov', filled circles with error bars) for the four bandpasses. Each color represents a different off-axis angle 
interval as annotated on the plot. Our best-fit models, as described in the Appendix, have been plotted as curves (see annotation); these models have been used to 
correct for incompleteness in our number count computations in § 2. 
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